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A COMPARATIVE STUDY OF THE DIGESTION OF PROTEINS 
AND CARBOHYDRATES IN GOATS DURING INFUSORIA- 
FREE AND -INFECTED PERIODS 


By Every R. Becxer,! J. A. ScHuLz,? AND M. A. EMMERSON? 
Iowa STATE COLLEGE 


Communicated August 6, 1929 


The purpose of the experiment was to determine the effect of the pro- 
tozoa upon the digestion of proteins and carbohydrates in ruminants. 
The goat was chosen for the experimental animal because it seems to 
lend itself to laboratory conditions better than our other common rumi- 
nants. ‘The experiment was planned so that a comparison of the results 
obtained in goats which had been rendered free of their rumen infusora 
for a definite period of time might be made with the results obtained for 
a like period after the animals were reinfected with the infusoria, which 
in our work were mostly members of the genus Entodinium and Diplo- 
dinium. The same amounts by weight of alfalfa hay and grain mixture 
were fed to any one goat during the infusoria-free period as during the 
period during which the goat carried the infusoria infection. Aliquot 
samples of the hay and grain which we fed to the animals were subjected 
to chemical analyses. By means of specially constructed crates all feces 
and urine from each goat were carefully collected. These materials were 
likewise subjected to chemical analyses. In all, four animals were con- 
ducted through digestion trials, and the digestibility coefficients of the 
protein, ether extract, nitrogen-free extract, crude fiber, hemicellulose, 
pentosans, alpha cellulose, and total dry.matter were determined. The 
partition of nitrogen between feces and urine, as well as the amounts of 
nitrogen stored daily, were also calculated. The results are shown in 
the accompanying table. 

An inspection of the above table will show that no significant advantage 
accrues to the host from the possession of the infusoria, so far as can be 
learned from the determinations made. This conclusion is quite con- 
trary to beliefs expressed by certain protozodlogists and physiologists. 
The complete results of the experiment are now being written up for 
publication at a later date. 
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THE MELANOPHORE ACTIVATOR OF THE EYE 


By BENJAMIN KRopp 
ZOGLOGICAL LABORATORY, HARVARD UNIVERSITY 


Communicated August 1, 1929 


The chromatic response of the skin to variation in background color 
depends mainly upon the presence of the eyes. In certain Crustacea 
(Koller, ’28; Perkins, ’28; Pouchet, ’72) and in some reptiles and am- 
phibians (Carlton, ’03; Kropp, ’27; Redfield, ’18) the eyes are the chief 
organs of light reception and as such mediate most light stimuli that 
ultimately affect the chromatophores. The fact that a blinded animal 
may change color slightly when transferred from a light to a dark back- 
ground does not alter the importance of the eyes, since in such cases the 
change from the expanded to the contracted condition, or in the reverse 
direction, is accomplished in twice the normal time and can be shown to 
be due to the presence of light receptors in the skin (Kropp, ’27; Laurens, 
"15; Parker, 03). In most cases, when an animal is totally blinded, the 
chromatophores assume a state of extreme expansion or contraction and 
at most show only slight variations from this condition. 

The nature of the réle that the eyes play in normal “adaptive’’ color 
changes is a moot question. Conceivably, they may function (1) as organs 
that transmit directly along nervous paths stimuli that affect the chroma- 
tophores, or (2) as organs that affect the chromatophores reflexly by in- 
ducing the secretion of chromatophore activating substances of the endo- 
crine glands. Perkins (’28) succeeded with Palaemonetes in obtaining 
the typical white-adapted condition in a blinded animal following the 
injection of a sea-water extract of eyes removed from white-adapted in- 
dividuals. ‘This was later confirmed by Koller (’28) on Crangon. It 
would seem, therefore, that there is produced in the eyes of these animals, 
under particular light conditions, a pigment-contracting substance which, 
carried in the blood stream, is responsible for the so-called white-adapted 
state. 

In the course of some work on the pigmentary responses of frogs and 
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tadpoles, an attempt was made by blood transfusion methods to test the 
hypothesis that in normal color changes the blood carried a melanophore 
activating material. While this had been successfully done with Phryno- 
soma (Redfield, 18), my experiments were negative. Later, the eyes of 
dark-adapted tadpoles of Rana clamitans were removed and extracted 
in distilled water. The extract was prepared and injected in the absence 
of light. When light was used the results were indifferent or negative. 
Such extracts were injected in small quantity into the abdominal cavities 
of normal light-adapted tadpoles and it was found that animals so treated 
showed the melanophore expanded state characteristic of dark-adapted 
tadpoles. This change occurred in from 5 to 30 minutes and lasted for 
2 to 3 hours regardless of light and background conditions. Four eyes 
were macerated quickly and thoroughly in 1 cc. of distilled water, filtered 
rapidly and 0.1 cc. to 0.2 cc. of the clear filtrate injected into the abdominal 
cavity. Controls were injected with an equivalent amount of distilled 
water. 

The reverse experiment, namely, the injection of extracts from the eyes 
of light-adapted animals into normal dark-adapted ones, has thus far 
yielded no clear-cut results. Nor has it been possible to obtain similar 
reactions in the case of adult frogs. However, with the minnow, Fundulus, 
the results have paralleled those with the tadpole and been even more 
satisfactory. Neither sea-water nor Ringer’s solution was of service as 
an extraction medium since they apparently destroyed in part the active 
material. Furthermore, the use of these fluids tends to complicate matters 
since each causes melanophore contraction when injected alone. 

The origin and nature of this melanophore activating agent cannot be 
stated with certainty. Preliminary work shows that the probable source 
is in the retina, since the reaction can be obtained when lens, cornea and 
iris are eliminated and the posterior half of the eye-ball alone extracted. 

Carlton, F. C., Proc. Amer. Acad. Arts and Sci., 39, 257-276, 1903. 
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NOTE ON MR. KING’S PAPER: “THE CRYSTAL STRUCTURE 
OF STRONTIUM” 


By F. Smwon anv E. VouSsEN { 
DEPARTMENT OF CHEMISTRY, SYRACUSE UNIVERSITY 


Communicated July 15, 1929 


In No. 4 of this periodical King stated that on redetermining the crystal 
structure of strontium he found a cubic face-centered lattice in agreement 
with our measurements,! but that the unit cubic edge (ao = 6.075 + 0.004 
A.U.) was 0.7% greater than ours, which figured 6.03 A.U. King believes 
that this rather considerable discrepancy is due to the fact that his measure- 
ments, which were carried out with a view to obtaining very accurate L 
results, were more exact than ours. We feel bound to point out that this { 
is not the reason, but that our spacing applies to liquid air temperature— 
a fact that King has apparently overlooked—whereas his value holds for 
room temperature. The difference in question cannot be calculated A 
accurately but can be gauged at 0.03 A.U. + 10% by the usual approxi- 
mation as shown in our paper. Computed for room temperature our 
result would thus be 6.06 A.U. The difference between this value and 
King’s is within the limits of our experimental error. We will not discuss 
the question whether our value or King’s is the more accurate, but we may 
mention that in a recent paper Ebert? gives the unit cubic edge of stron- 
tium as 6.05 A.U. at room temperature. 


1 Zeit. phys. Chemie, 133, 165, 1928. 
2 Zeit. anorg. Chemie, 179, 4, 1929. 


THE RAMAN SPECTRUM OF GYPSUM 


By Roscog G. DICKINSON AND RoBERT T. DILLON 
GaTtEs CHEMICAL LABORATORY, CALIFORNIA INSTITUTE OF TECHNOLOGY 


Communicated July 27, 1929 


In an investigation of the Raman spectra of some ionized substances 
in aqueous solution, it was found! that the frequency shift given by a 
solution of potassium carbonate was in fairly close agreement with the 
strongest shift given by calcite. Data for such a comparison in other 
cases do not appear to exist, for change of frequency on scattering has 
been measured for only a few crystalline substances. In addition to 
calcite,?*> these are quartz,?***® ice,> and topaz.’ No modified fre- 
quencies were found in the light scattered by halite,? NaCl, fluorite,® 
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CaF:, or by the salts’ LiF and NaF. In the present paper are described 
the results of measurements on gypsum, CaSO,.2H.O, and a comparison 
of these with data on ammonium sulfate solution. 

The crystal used was a clear, colorless specimen internally free from 
cleavage cracks. It was about 15 cm. long and had a cross-section of 
about 2 cm.?; the long edges of the prism faces were parallel to the c-axis. 
The crystal was placed in a water-jacketed glass tube with an end window 
and irradiated with light from a mercury arc. The arc was placed parallel 


TABLE 1 
RAMAN SPECTRA FROM GYPSUM 
FREQUENCY OF INTENSITY OF WAVE LENGTH OF FREQUENCY OF FREQUENCY 

SHIFTED LINE, SHIFTED LINE EXCITING LINE, EXCITING LINE, CHANGE, 

cm. ~! cm. ~! cm, ~! 
24291 .5 Weak 4046 .56 24705 .5 414.0 
24085 .0 Very weak (broad) 4046 .56 24705 .5 620.5 
24033 .6 Very weak 4046 .56 24705 .5 671.9 
23985 .7 Medium (broad) 3650.15 27388 .4 3402.7 
23949 .3 Weak (broad) 3654 .83 27353 .3 3404.0 
23896 .1 Strong (very broad) 3650.15 27388 .4 3492.3 
(3663 .27) (27290 .5) (3394 .4) 

23361 .4 Medium weak (broad) 3654 .83 27353 .3 3491.9 
23798 .1 Medium weak (broad) 3663 .27 27290 .5 3492 .4 
23697 .8 Strong 4046 .56 24705.5 1007 .7 
23569 .3 Medium 4046 .56 24705 .5 1136 .2 
23507 .2 Very weak 4077 .8 24516.1 1008 .9 
23351 .5 Very faint 4358 .34 22938 . 1 —413.4 
22523 .0 Medium weak 4358 .34 22938 .1 415.1 
22445 .2 Very weak 4358 .34 22938 .1 492.9 
22318 .6 Weak (broad) 4358 .34 22938 .1 619.5 
22267 .3 Weak 4358 .34 22938 .1 670.8 
21987 .1 Very weak 4347 .50 22995 .3 1008.2 
21929 .9 Very strong 4358 .34 22938 .1 1008 .2 
21803 .3 Medium strong 4358 .34 22938 .1 1134.8 
21303 .2 Medium strong (broad) 4046 .56 24705 .5 3402 .3 
21215.2 Strong (broad) 4046 .56 24705 .5 3490.3 


to the c-axis of the crystal, and scattered light from the crystal was taken 
off in the direction of the c-axis. Since the end of the crystal was quite 
irregular, the crystal was submerged in chlorobenzene which has a refrac- 
tive index close to those of gypsum. ‘To insure that no scattered lines 
had reached the spectrograph from the liquid, carbon tetrachloride, 
which gives a quite different Raman spectrum from chlorobenzene, was 
_ used in one case with no change in the observed scattered lines. The 
glass tube containing the crystal was painted black around the ends of 
the crystal. The light from the crystal was converged with lenses and 
thrown on the spectrograph slit with the aid of a reflecting prism. 

The spectrograph was a two-prism glass instrument arranged with 
interchangeable camera lenses so as to give either a dispersion of about 
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100 cm.~! per mm. with an aperture of f 8, or a smaller dispersion and 
higher speed corresponding to an aperture of f 3. A comparison spectrum 
of the iron arc was recorded on each plate, and the plates measured on a 
Gaertner comparator. The frequencies of the shifted lines were computed 
by linear interpolation between the frequencies of closely neighboring 
iron lines. 

The measured frequencies and estimated intensities of the modified 
lines are given in the first two columns of table 1. ‘The wave-lengths and 
frequencies of the corresponding exciting lines are given in the third and 
fourth columns, respectively, while the frequency shifts are given in the 
last column. As a control on the assignment of exciting lines to the 
various observed lines, we have made a photograph while flowing through 
the water-jacket an alkaline solution of potassium chromate of a con- 
centration such that considerable light of wave-length \ 4358 was trans- 
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mitted, while that of shorter wave-length was largely absorbed. All of 
the modified lines appearing on this photograph were ones for which 
excitation by \4358 is indicated in table 1, and all lines thus indicated 
appeared on the photograph with the exception of the very faint anti- 
Stokes line. 

The data of table 1 indicate eight different frequency shifts; the best 
values of these are summarized in table 2. 

TABLE 2 
FREQUENCY SHIFTS GIVEN BY GYPSUM 


Frequency 
shift (em.~!) 414.2 492.9 620.0 670.4 1008.2 1185.5 3403.0 3491.5 
Intensity Med. Very Weak Weak Very Med. Med. Strong 
weak weak (broad) strong strong strong (broad) 
(broad) 


We have also rephotographed the spectrum from ammonium sulfate 
solution using higher dispersion than previously.'_ The new measurements 
yield the following values for the frequency shifts of the Raman lines pre- 
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viously reported: 451, medium (broad); 620, medium (broad); 980.3, 
strong. In addition there appears to be a very broad shift with Av ~ 
1113 as well as the still broader band due to water with Av equal to about 
3430. With this higher dispersion it is especially evident that the strong- 
est line is far sharper than the others. A comparison of this spectrum 
with that of gypsum is shown in figure 1. 

The two largest frequency shifts given by gypsum are presumably due 
to water of crystallization while the smaller shifts are very likely due 
to sulfate ion. Schaefer and Schubert* measured residual rays from a 
variety of sulfates and concluded that all sulfates possess prominent 
maxima of reflection in the neighborhood of 9 and 16u which are to be 
attributed to the sulfate ion. In the case of gypsum, the maxima occurred 
at 8.74 and 14.84 with a weak maximum at 8.62u, and another weak 
one at about 3u due to water of crystallization. The two prominent 
maxima correspond to frequencies of 1143 and 674 cm.~', respectively, 
in reasonable agreement’ with the values 1135.5 and 670.4 which we have 
found. The two smallest Raman shifts do not fall in the range of the 
residual ray measurements. The absence of any maximum corresponding 
to the strongest Raman line (Av = 1008.2) suggests that the initial and 
final states of the transition producing it are without electric moment, 
as Carelli, Pringsheim and Rosen" have suggested for an analogous line 
in the case of calcite. This appears to afford a reasonable explanation 
of the sharpness of this line in aqueous sulfate solution; for absence of 
electric moment would presumably minimize the interaction between 
sulfate ion and water molecules. 

If the symmetry of the sulfate ion in aqueous solution may be taken as 
tetrahedral, Brester’s'' considerations concerning the multiplication of 
characteristic frequencies consequent upon removal of symmetry afford 
one possible explanation of the multiplicity of lines found from gypsum. 
The point group symmetry of this crystal is monoclinic holohedral. ‘The 
space groups symmetry has been reported by Onorato'® to be C3, with 
8CaSO,.2H2O in a unit of structure that employs a face-centered lattice; 
if an end-centered lattice is used, the number of CaSO;.2H.,O becomes 
four. If the sulfur atoms occupy equivalent positions, the sulfate ion can, 
strictly speaking, have only either a center of symmetry, a two-fold axis, 
or a plane of symmetry in this crystal although it may approximate to a 
higher symmetry. If the sulfur atoms are not all in equivalent positions, 
a multiplication of lines could result for this reason. In any case, it is 
also possible that the weaker Raman lines do not represent fundamental 
frequencies. 


1 Dickinson and Dillon, Proc. Nat. Acad. Sci., 15, 334 (1929). 
2R. W. Wood, Phil. Mag., 6, 729 (1928). 
3 Gr. Landsberg and L. Mandelstam, Zeit. Phys., 50, 769 (1928). 
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4 P. Pringsheim and B. Rosen, Zeit. Phys., 50, 741 (1928). 

51. R. Rao, Indian J. Physics, 3, 123 (1928). 

6H. Nisi, Proc. Imp. Acad. of Japan, 5, 127 (1929). 

7C. Schaefer, Zeit. Phys., 54, 153 (1929). 

8 Schaefer and Schubert, Ann. Phys., 50, 283 (1916). 

9 This agreement may be fortuitous because of the fact that maxima of reflection 
and absorption are not quite coincident and are sometimes very appreciably removed 
from each other. See, for example, Havelock, Proc. Ray. Soc. 105A, 488 (1924). 

10 Carelli, Pringsheim and Rosen, Zeit. Phys., 51, 511 (1928). 

11 Brester, Jbid., 24, 324 (1924). 

12, Onorato, Zeit. Krist., 66, 504 (1928). 


RAMAN SPECTRA FROM ACETONE 


By Rosert T. DILLON AND Roscog G. DicKINSON 
GaTEs CHEMICAL LABORATORY, CALIFORNIA INSTITUTE OF TECHNOLOGY 


Communicated August 1, 1929 


Measurements of Raman spectra from acetone have been published in a 
recent paper.! These measurements were so interpreted that out of 
thirteen modified lines, twelve corresponded to frequencies of observed 
infra-red absorptions, the thirteenth lying outside the range of the infra- 
red measurements. This result is surprising; for the occurrence of fre- 
quency shifts in Raman effect without the appearance of the corresponding 
infra-red absorptions is a very common phenomenon. We have examined 
Raman spectra from acetone; and as both our measurements and our 
interpretation disagree with those mentioned above, it seems desirable to 
report them. 

The acetone used was Merck’s ‘‘C.P.”’ It was distilled through a 30 cm.- 
Eastman fractionating column the product used all coming over within 
0.05° of 55.7° at 746 mm. It is important to use materials of reasonable 
purity; we have, for example, found it possible to bring out the strongest 
benzene lines from a carbon tetrachloride solution containing only 0.46% 
benzene. 

The acetone was irradiated with light from a glass mercury arc, and the 
spectra photographed with a glass spectrograph having a dispersion of about 
100 cm.~! per mm. in the blue. The iron arc was used as a standard. 

The measured frequencies of the shifted lines together with their in- 
tensities are given in the second and third columns of table 1. Com- 
parison of these measurements with those which have already been re- 
ported! shows considerable disagreement even among the strongest lines. 
Apparently seven lines are concerned in the production of these Raman 
lines. The wave-lengths and frequencies of these mercury lines are: 
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I 3650.15 A 27388.4 cm.-! V 


II 3654.83 A 27353.3 em.-? 
III 3663.27 A 27290.5 cm.-! 


IV 4046.56 A 24705.5 cm.-! 


TABLE 1 


DILLON AND DICKINSON 





Proc. N. A. S. 


4077.8 A 24516.1 cm.~ 


OBSERVED RAMAN LINES FROM ACETONE 


OBSERVED 


FREQUENCY OF 
SHIFTED LINE, 


cM. ~ 


24464. 
24428. 


24383 . 


24367 . 


24316. 


24284. 
24212. 
24174. 
23918. 
23727 
23639 
23482. 
23366 . 
23275. 
22546. 
22446. 
22407. 
22207 . 
22150 
22032. 
22009 
21869. 
21782 
21738. 
21712 
21696 . 
21590. 


21507. 


21227 
21190. 
20069 . 
20015 
19972. 
19931. 


OQNaAnMANI SCH Ww 


moO MO 


1 
1 
9 
oO 
4 
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4d 


a 
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Na © 


> 
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ESTIMATED 
INTENSITY OF 
SHIFTED LINE 


Very strong 
Strong 


Medium 
Strong 
Weak 


Weak 

Very weak 
Medium weak 
Strong 
Medium weak 
Medium weak 
Medium 

Very weak 


Medium (broad) 


Medium 

Weak 

Medium 

Very weak 
Very strong 
Very weak 
Very weak 
Medium 

Very strong 
Medium strong 
Medium strong 
Medium strong 
Medium strong 


Medium strong 
(broad) 
Strong 
Very weak 
Very weak 
Strong — 
Medium weak 
Medium weak 


ASSUMED 
ORIGIN OF 
LINE 
I—m 
I-—n 
Il — m 
I-o 
II-—n 
III — m 
II —o 
III —n 
IV -—a 
III — o 
IV — dD 
IV -—c 
Iv -—d 
V-d 
IV —g 
IV -—h 
IV -—14 
IV -j 
VII —a 
VII — b 
VII —c 
VI -—d 
VII -—-d 
VII —e 
VII —f 
VII — g 
IV —m 
IV —n 
VII -—h 
IV -—o 
V—m 
VII -—i 
VII — j 
VII —k 
VII —/ 
VI-—m 
VII — m 
VII —n 
VII —o 


VI 4347.50 A 22995.3 cm.-! 
VII 4358.34 A 22938.1 cm.-! 


CALCULATED 
FREQUENCY OF 
SHIFTED LINE, 
cm. ~! 


24464. 
24421 
24429 . 
24382 . 
24386 . 
24366 . 
24346. 
24323 . 
24314. 
24284 . 
24213. 
24175. 
23917 . 
23728 
23638 . 
23482. 
(23366 . 
23275. 
(22546. 
22446 . 
22408 .1 
22207 .3 
22150.1 
(22032 . 5) 
(22009 .7 
21870. 
21781. 
21738 .9 
21715.3 
21699 .1 
21592 .3 
21599 .1 
21507 .6 
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Ni 


(21227 .1) 
(21190. 1) 
20071 .5 
20014.3 
19971.5 
19931 .7 
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The strongest of these lines are I, IV and VII. Lines V and VI are so 
weak as to be involved in the production of only the most prominent 
Raman shifts. The modified lines are all accounted for by using the 
frequency differences summarized in table 2. 


TABLE 2 
FREQUENCY SHIFTS FROM ACETONE 
FREQUENCY RELATIVE RELATIVE 
DESIGNATION SHIFT, INTENSITY OF \ RAMAN ALR. INTENSITY OF 
OF SHIFT cm.~! SCATTERING BM ABSORPTION 
a 391.5 Medium (25.6) 
b 492.0 Weak (20.3) 
c 530.0 Medium (18.9) 
13.10 2.3 
d 788 .0 Very strong 12.69 12.63 2.2 
e 905.6 (?) Very weak 11.08 11.30 12 
f 928.4 (?) Very weak 10.77 10.90. 4.0 
10.40 0.2 
9.70 0.2 
g 1067 .4 Medium 9.38 9.10 2.0 
h 1222.8 Medium strong 8.17 8.30 4.0 
1 1339.0 (?) Very weak 7.47 7.60 0.3 
7.30 0.5 
j 1430.5 Medium (broad) 6.99 6.98 0.4 
6.60 0.5 
k 1711.0 Strong 5.84 
l 1748.0 (?) Very weak 5.72 5.75 5.2 
4.70 1.5 
m 2923 .8 Very strong 3.42 3.41 4.0 
n 2966 .6 Strong 3.37 
0 3006 .4 Strong 3.27 


In the fifth column of table 1 are given calculated frequencies of modified 
lines obtained assuming the various lines to originate in the way indicated 
in the fourth column. (For example, I—m is the mercury line I of fre- 
quency 27388.4 cm.~! shifted by the amount 2923.8 cm.~! which is shift 
m.) It may be seen that the assumed origins of the lines are such that 
modified lines corresponding to any one shift have qualitatively (allowing 
for the smaller sensitivity of the plate toward the green) the same relative 
intensities as the mercury lines producing them. Several of the lines in 
table 1 (Nos. 2, 3, 4, 5 and 27) should result from the approximate super- 
position of two or more modified lines; measurements of these have been 
disregarded in selecting the best values given in table 2. 

The fairly strong shift k (Av = 1711.0) occurs only once (line 29); 
Hg 4046 shifted by this amount would produce a line almost coincident 
‘with Hg 4347. If, however, a different origin were assumed for line 29 
we should expect to find other lines not occurring on the plate. To 
account for lines 13, 15, 20, 21 and 30 corresponding to shifts 7, a, e, f and 
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l we have also found it necessary to assume five frequency shifts not other- 
wise employed. When a potassium chromate filter was used, line 15 
was greatly strengthened compared with, for example, lines 12 and 14; 
this leaves little doubt that line 15 arises from Hg 4358 modified by shift 
a. There remains, possibly, some ambiguity in the origin of lines e, f, z 
and/. We have attempted to remove this ambiguity by using the potas- 
sium chromate filter; but because of the extreme weakness of these lines 
we have not been entirely successful. 

In the fourth column of table 2 are given calculated infra-red wave- 
lengths corresponding to the frequency shifts. In the fifth and sixth 
columns are the wave-lengths and intensities of the maxima of infra-red 
absorption observed by Coblentz.? If a latitude of 0.34 be allowed, all 
of the observed Raman shifts can be brought into agreement with infra- 
red absorption maxima in so far as the investigated regions overlap. In 
view of the considerable number of frequencies involved in both spectra, 
the necessity of allowing so much latitude leaves some of the coincidences 
unconvincing. In any case, the intensities in scattering and in absorp- 
tion bear no evident relation. This question has been discussed elsewhere.’ 

1J. W. Williams and A. Hollaender, Proc. Nat. Acad. Sci., 15, 421 (1929). 


2 Coblentz, “Investigation of Infra-red Spectra,’ Carnegie Inst. Publications, 1 (1905). 
3 Dickinson, Dillon and Rasetti, Phys. Rev., 34 (1929). 


THE THERMAL DECOMPOSITION OF ACETONE IN THE 
GASEOUS STATE 


By F. O. Rick anv R. E. VOLLRATH* 
DEPARTMENT OF CHEMISTRY, THE JoHNS HopKINs UNIVERSITY 


Communicated August 2, 1929 


In the field of reaction kinetics, the study of unimolecular homogeneous 
gas reactions is, from a theoretical standpoint, most interesting. At the 
present time but few really homogeneous, unimolecular reactions are 
known and this work was undertaken with the object of extending the 
number of organic decompositions of this type. In all the investigations 
of such decompositions made so far, the method used to follow the reaction 
has been to measure the rate of increase of pressure and from this to ob- 
tain the velocity constant. This method has two obvious disadvantages: 
(1) we must know what the products of decomposition are, and (2) these 
must not decompose or react with each other appreciably at the tempera- ° 
ture of the experiment. Unfortunately in the case of organic compounds, 
these conditions seriously limit experimental work so that the decomposi- 
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tion of an organic series cannot be followed using this method. Conse- 
quently we have developed a method in which we follow the concentra- 
tion of the decomposing compound analytically and in this way avoid 
all the troublesome uncertainties connected with the products of de- 
compositions; provided that these do not react with the decomposing 
substance, and that we are dealing with.a single homogeneous unimolecular 
change, our measurement of the rate should be free from error. ‘The 
limitation of this method is that we must devise a quantitative method for 
estimating the organic compound; fortunately this can be done in a great 
many classes of compounds so that this method enables us greatly to 
extend the range of these investigations. 

1. Previous Work.—Hinshelwood and Hutchinson' have recently mea- 
sured by a manometric method, the rate of decomposition of acetone vapor 
in a quartz bulb maintained at temperatures of 500-600°C. They de- 
duced from their results that the decomposition is unimolecular and in 
accordance with the following scheme: 


CH;COCH; —> CO + CH; + CHs. 


The two methyl groups then react to give methane, ethylene and hydrogen. 

It is well known’, however, that when acetone is passed through a heated 
glass tube, methane separates from the molecule and ketene results, ac- 
cording to the following equation, 


CH;COCH; —> CHe:CO + CH, 


and Taylor*® has suggested the possibility that Hinshelwood and Hutchin- 
son really measured the rate of thermal decomposition of ketene. Apart 
from any uncertainty in the previous work it seemed desirable to develop 
a method which would be applicable not only to acetone but to its homologs, 
since a knowledge of their rates of decomposition and energies of activation 
would be of great interest. 

2. Experimental Method.—A flow system was employed in which pure 
dry nitrogen saturated with the vapor of the acetone or higher ketone 
was passed at a known rate through a quartz tube heated in an electric 
furnace. The gases were then passed into a series of small wash bottles 
containing water in which the acetone could be estimated by hydroxyl- 
amine hydrochloride,‘ or by Messinger’s method.® This method had 
been already worked out in this Laboratory;* in order to prove that the 
products of the decomposition did not react with the acetone it seemed 
desirable to show that one molecule of CO was present for every molecule 
of acetone decomposed. Dr. Slagle’ had shown that the quantity of CO 
produced was much less than that expected. Consequently we tested 
for the presence of ketene in the gases leaving the tube and found very 
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considerable quantities. We therefore made measurements of the amounts 
of ketene formed when acetone is heated in a quartz tube. 

3. Experimental Results ——A preliminary experiment confirmed the 
presence of ketene in the gases coming from the quartz tube at 630°C.; 
these were passed through a solution of aniline in ether and the excess 
aniline and ether was evaporated off. The.residue was shown to be 
acetanilide by the melting point and by mixing the residue with acetanilide 
and taking a second melting point. In this preliminary experiment about 
20% of the acetone introduced into the quartz tube was converted into 
ketene. The amount of ketene was also determined by absorbing it in 
N/10 sodium hydroxide and measuring the amount of base neutralized. 
Using this method, a number of runs were made to determine the influence 
of rate of flow on the yield of ketene and the results are summarized in 
table 1. An additional series of experiments was then performed to find 


TABLE 1 
‘TEMPERATURE OF FuRNACE = 635°C. 
Min. in furnace 1.56 1.56 2.36 2.36 4:00 4.00 4:00 


Mols ketene formed per 
mol. acetone introduced 0.256 0.260 0.255 0.254 0.249 0.254 0.249 


what proportion of the decomposed acetone could be recovered as ketene 
and the results are given in table 2. The average half-life for acetone at 


TABLE 2 
TEMPERATURE OF FURNACE = 588°C. 
Acetone introduced (gram) 0.4135 0.4213 0.3357 0.4008 
Acetone undecomposed (gram) 0.3522 0.3624 0.2711 0.3499 
Acetone decomposed (gram) 0.0613 0.0589 0.0646 0.0509 
Mols ketene formed per mol. 
acetone (decomposed) 0.559 0.582 0.552 0.622 

Time in furnace (min.) 2.42 2.45 2.47 2.42 


this temperature is 2.5 min. as compared with 3.1 min. obtained by 
Hinshelwood and Hutchinson.! ; 

4, Discussion.—The results obtained so far show that for every 100 
molecules of acetone decomposed approximately 60 molecules of ketene 
can be recovered, which is about the same as the yields reported by Hurd 
and Tallyn’ who heated the acetone ina pyrextube. It seems likely, there- 
fore, that the primary process in the homogeneous unimolecular decom- 
position of acetone is the separation of a molecule of methane and forma- 
tion of ketene. Presumably at the high temperature in the furnace the 
ketene undergoes a bimolecular decomposition as follows: 


2CHe:CO = CoH, + 2CO. 


According to this scheme 100 molecules of decomposed acetone should 
give 100 molecules of methane, 30 molecules of ethylene and 60 molecules 
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of carbon monoxide, which corresponds very closely with the results of 
the gas analysis made by Hinshelwood and Hutchinson.' 

5. Conclusion—The decomposition of acetone has been followed by 
an analytical method and it has been found that about 60% of the acetone 
decomposed can be recovered as ketene. Methods in which rate of change 
of pressure are used for following the course of the reaction are therefore 
unreliable. 


* Du Pont FELLOW. 

1 Hinshelwood and Hutchinson, Proc. Roy. Soc., 111A, 245 (1926). 

2 See Hurd, The Pyrolysis of Carbon Compounds, p. 247. Chemical Catalog Co., 
1929. 

3 Taylor, J. Phys. Chem., 30, 1433 (1926); book review. 

4 Meyer, Lehrbuch der Organisch-Chemischen Methodik, p. 846. Springer, 1922. 

5 Goodwin, J. Am. Chem. Soc., 42, 39 (1920). 

6 Slagle, Dissertation, ‘‘Decomposition of Aliphatic Ketones in the Gas Phase,” 
Johns Hopkins University, 1929. 

7 Hurd and Tallyn, J. Am. Chem. Soc., 47, 1427 (1925). 


THE EQUILIBRIUM CONSTANTS OF REACTIONS INVOLVING 
HYDROXYL 


By DONALD STATLER VILLARS 
CHEMICAL LABORATORY, UNIVERSITY OF ILLINOIS 


Communicated July 26, 1929 


A recent interpolation of band spectra data has indicated that the 
heat of dissociation of OH should be 138 k. cal.' An accurate value for 
the heat required in this dissociation would enable one to calculate the 
equilibrium constants of reactions involving this molecule.?, Thus Bon- 
hoeffer® and Reichardt have calculated the equilibrium constant of the 
following reaction, using their value for the heat of dissociation of water 
into H and OH: 

_ [H] (OH? 


2H.2O = He + 20H Kou = 
2! 2 OH [HO]? 


As my value for the heat of dissociation of OH is some 14 k. cal. larger 
than that corresponding to Bonhoeffer and Reichardt’s data, I have re- 
calculated their constants and contributed some additional ones. 

In order to calculate the free energies of these reactions, use was made 
of the familiar equations: 


F=H-—TS 


AF = AH — TAS = —RT Ink 
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AH = Do os > (3 RT + Evot. + vibr. +, 9 


AF _ Do , 
: beitaille deakdienas 


where 
— Do = Qo = heat of reaction at absolute zero and 


ae ee 4 ©. 
S=S'+ 5) + Tr 

The entropies were obtained from an application of the method of 
Gibson and Heitler.t In these methods involving the statistical inter- 
pretation of the third law of thermodynamics, the entropy of the mole- 
cule is considered to be the sum of its translational entropy, 


9 S/ety é 
S = Rig Cee = 

eee 
its rotational entropy, its vibrational entropy, and that due to all its 
different possible kinds of orientation. The latter may be expressed® in 
short by 


Sm = = + Rin > p,e~ "7, 


where /, is the a priori probability of the state n. This applies equally 
well to entropy of vibration and rotation, and, in these cases, E,, is the 
total amount of the corresponding energy, i.e., Epitration Evotation orientation 
(= 0 for a degenerate state). 

As a further check on the reliability of this method, beyond that already 
shown by Gibson and Heitler, the entropies of oxygen and hydrogen at 
25°C. were calculated. The count of the rotational entropy is simplified 
in the case of oxygen because there is only one kind present (symmetric 
or antisymmetric), as is shown by the fact that its band spectrum con- 
sists of alternately missing lines. Its theoretical entropy comes out to 
be 50.7, whereas the thermal values range from 45.6 to 49.2.6 This high 
value comes from the inclusion of the *S spectroscopic state of the oxygen 
molecule.’ If it were not for its triplet character, the calculated entropy 
would be 48.5. If the discrepancy is not due to calorimetric errors, this 
value predicts a fourth allotropic form of solid oxygen to be found with 
a transition temperature below 17°K. 

The count of the rotational entropy of hydrogen is complicated by the 
presence of the symmetric and antisymmetric kinds in a one to three 
ratio. The question has been raised by Giauque* whether the entropy 
of mixing of these two kinds should not be included, especially as the 
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two forms do not intercombine spectroscopically. Rodebush*® has an- 
swered that there really is only one hydrogen and that although the two 
forms do not intercombine in their spectral relations, they do so by col- 
lision when given a sufficient length of time to reach equilibrium. It is 
comforting to know that it really makes no difference in the final result 
which method of counting is used. The entropy of hydrogen at 25°C. 
comes out to be 34.0'°, whether we consider it two kinds and include the 
entropy of mixing or one substance which will reach equilibrium if given 
time enough. This, in fact, must be,'! as can be seen from the following: 


40? & kT 7 


par = Kin (3 + 1) + 








(one substance Gibson and Heitler count) or 


Sree ad Santi + Sign + Sentetes 











= nat ml n3™ oe 1 1 
3 3 Exot 
4m? TkT Fue 
= Rit me 4+ ak 


The entropy of hydrogen was also calculated by a summation process 
to check against possible inaccuracies introduced by the integration of the 
rotational sum. ‘The results were identical. 

The entropy of OH was next computed. As this molecule does not give 
rise to nuclear degeneracy the number of rotational terms is not halved as 
in the case of hydrogen and oxygen and the rotational entropy is expressed 
by the formula: 

Sro = Rin 8x? IRT/h? + E,.,/T. 


In addition to this rotational entropy comes entropy to account for the 
electronic states in the *P inverted ground term, which was incorrectly 
accounted for by Bonhoeffer and Reichardt. This should be as follows: 


Sate = Kin (4 + he) + ae 


The upper doublet term is not negligible in the case of the OH molecule, 
nor is it of full weight. 

Table 1 gives the numerical data used to calculate the entropy formulas 
given on the following page. 
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TABLE 1 
I wo Av 
Hydrogen 4.67 X 10-*! g. cm.? 4415.2 cm.7! 
Oxygen 19.2 x 10-* 1565.4 
Hydroxyl 1.43 X 10-® 3570 137.9 cm.~! 
So, = 13.70 + 11.52 log T — 4.61 log (1 — e777) + 2 R 
: Eo, 
— 461 log Patm + - 
Su, = —1.51 + 11.52 log T — 4.61 log (1 — e~9/7) +2 R 


— 4.61 log pom + Ew 





5 

Sou = 5.77 + 11.52 log T — 4.61 log (1 — e~*!!”/7) 
, 9,—197.5/T a) i Fou 
+ 4.61 log (4 + 2e )+5R — 461 log parm + a 


From these values and the free energy formulas from Lewis and Randall, 


AF °i3,0 gas = —57410 + 0.94 T In T + 0.00165 T? 
—0.00000037 T* + 3.92 T 
(Lewis and Randall, 485) 


AF °v10, gas = — 31200 + 5.5 Tin T — 0.00115 T? — 9.3 T 
(Lewis and Randall, 496) 


the equilibrium constants of the following reactions were calculated: 





20H = He “+ Oz . (1) K, oe [He] [Oz] 
{OH}? 
, [He]? [Oz] 
2H = 2H. + Oo. — 114k. cal. 2) Ko, = —— 
sieeueee 
[He][OH]? 
2 HO = H. + 20H 3) K se 
@) Kou = Fle 
20H = HO, (4)K, = [202], 
[OH]? 


The results are given in table 2, in which the subscripts B denote that 
Bonhoeffer and Reichardt’s values of the heat of dissociation of water were 
used and the subscripts V denote that my value was used. 
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TABLE 2 
TEMPER- LOG Ko, LOG Ko, LOG Koy 
ATURE LanpR, BanvdR, B anv R, 
ABS. p. 485 p. 96 LOG Kon, p. 96 LOG Kony L0G Ki, Loc Kiy Loc Kug Loc K,y 


1000 —-19.8 —20.1 —21.8 —21.1 —15.5 +1.95 —4.35 +2.74 —3.55 

1300 —-13.9 —14.01 —15.2 —14.3 —10.34 +1.31 —3.53 +0.30 —4.54 

1500 —11.2 —11.42 —12.24 —11.4 — 8.05 +1.03 —3.16 —0.79 —4.98 

1705 — 9.13 — 9.28 — 9.95 — 9.0 — 6.26 +0.81 —2.87 —1.66 —5.35 
7.2 


1900 -— 7.55 — 7.6 — 8.21 — — 4.90 +0.65 —2.66 —2.29 —5.60 
interp. 

2155 — 5.94 — 6.08 — 6.42 — 5.5 — 3.50 +0.48 —2.43 —2.99 —5.91 

2505 — 4.27 — 4.31 — 4.59 — 3.5 — 2.08 +0.32 —2.19 —3.67 —6.18 


20H —> Hz + O, + 14,000 cal. (B) 
20H —_—> He + O, et 15,000 cal. (V) 

1 Villars, J. Am. Chem. Soc., 51, 2374-7 (1929). 

2 Dr. Langmuir has kindly pointed out that the knowledge of such equilibrium 
constants will be extremely useful in investigating the possibility of the direct formation 
of hydroxyl during the combustion of hydrogen and oxygen at high temperatures, an 
investigation which he is at present carrying on. 

3Z. physik. Chem., 139A, 75-97 (1928). 


4Z. Physik, 49, 465. 
5 Cf. Rodebush’s chapter XVII, p. 1202, in Taylor, A Treatise on Physical Chemistry, 


D. Van Nostrand, 1925. 
6 Lewis and Gibson, J. Am. Chem. Soc., 39, 2554 (1917). 
7 Mulliken, Phys. Rev., 32, 186. 
8 J. Am. Chem. Soc., 50, 3221. 
§ These PROCEEDINGS, 15, 678-680 (1929). 
10 Rodebush, loc. cit., obtained this value and showed that the experimental method 


used to determine the entropy of hydrogen is theoretically unsound. 
11 Prof. Rodebush has shown independently that this relation must obtain. 


ON THE CRYSTAL STRUCTURE OF THE CHLORIDES 
OF CERTAIN BIVALENT ELEMENTS 


By Linus PAULING 
GaTEs CHEMICAL LABORATORY, CALIFORNIA INSTITUTE OF TECHNOLOGY 
Communicated August 6, 1929 


In 1925 Bruni and Ferrari! reported that the lines on a powder photo- 
graph of the hexagonal crystal MgCl, could be accounted for on the basis 
of a rhombohedral (pseudocubic) unit of structure with a = 5.08 A and 
a = 90° containing 2 MgCl. Later they announced? that ZnCh, CdCl, 
and MnCl, have the same structure, with values of a a few degrees larger 
than 90°. Ferrari? then added CoCl, and NiCl, to the group; certain 
weak lines occurring on the powder photographs of these substances re- 
quire the value of a to be doubled, giving a rhombohedral unit with a = 90° 
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and a = 10.02 A and 10.00 A, respectively, containing 16 M+*+Cl. Gold- 
schmidt‘ has also reported that unpublished powder photographic data 
obtained by Oftedal in his laboratory show that RuCh, RhCh, PdCl:, 
IrCl, and PtCl, are to be assigned the same structure as NiClo. 

None of these investigators succeeded in deducing the atomic arrange- 
ment in these crystals; the utilization of Bruni and Ferrari’s data to this 
end is described in the following paragraphs. 

Bruni and Ferrari published reproductions? of Laue photographs of 
CdCh, MgCl, and MnCk, taken in order to determine the axial ratios for 
these crystals. I have assigned indices to the spots appearing on the best 
of these photographs, that of CdCl, with the aid of a gnomonic projection 
made from a tracing. The forms producing these spots and values of nd 
calculated for a unit with a = 94° and a = 10.2 A (an estimated value, 
about equal to that for ZnCl.) are given in the first and second columns 
of table 1. 

Only a few forms show m) values below 0.4 A, indicating first-order re- 
flections. It is, indeed, possible to account for the production of all of the 
spots on the photograph on the basis of a much smaller unit, a rhombo- 
hedron with a = 6.35 A and a = 36° 40’, containing 1MCk. This unit 
has lattice points at the face-centered positions of the rhombohedron 
which is itself formed by the face-centered positions of the original pseudo- 
cubic unit. Inasmuch as values of md calculated for this unit (fourth 
column of table 1) in no case fall below the short wave-length limit, the 
unit may be accepted as the true one. 

The presence of a three-fold axis and three planes of symmetry in the 
Laue photograph taken with the incident beam normal to the basal plane 
requires that the crystal be isomorphous with one of the point-groups 
Cy, D3 and D3z. The fundamental lattice has been found to be rhombo- 
hedral. The only ways of placing 1CdCl, in the unit compatible with 
these restrictions are 


I Cd at 000, 2 Cl at uuu, fai; Dj and D},; 
Il Cd at i, 2 Clat uuu, aaa; Dj and D3,; 
III Cd at uuu, Cl; at vvv, Cly at www; C},. 


In the absence of any evidence requiring it, the three-parameter structure 
III need not be considered. 

Structures I and II are identical. The data are not sufficient for an 
accurate determination of the value of the chlorine parameter u. The 
fact that the powder photograph of MgCl: is nearly identical with that of 
LiCl indicates that the chlorine atoms are in nearly the same positions 
in the two crystals, namely, at the points of a face-centered cubic lattice. 
This corresponds to a value of u of '/,. Further evidence that the chloride 
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ions are in approximate cubic close-packing is provided by the considera- 
tion of interionic distances. With a crystal radius of 1.81 A,* 32 chloride 
ions would occupy a cube 10.24 A on an edge; this approximates closely 
the values reported for the various bivalent chlorides. 





TABLE 1 
DATA FROM BRUNI AND FERRARI'S LAUE PHOTOGRAPH OF CADMIUM CRLORIDE 
PSEUDOCUBIC UNIT TRUE UNIT ESTIMATED STRUCTURE 
a= 102A a = 94° a = 6.36 A, a = 36°40’ INTENSITY* FACTOR 
{hkl} nd {hkl} ny 
335 0.33 A 111 0.33 A m Cd — 2C1 
331 0.68 211 0.34 w. Cd — 21 
557 0.40 221 0.40 m. Cd 
122 1.44 ° 321 0.36 m.s. Cd + 201 
779 0.32 331 0.32 v.w. Cd 
320 1.04 411 0.26 V.w. Cd — 201 
731 0:72 631 0.36 v.w. Cd — 2C1 
411 1.52 611 0.38 m. Cd + 2Ci 
432 140. 750 0.35 m.w. Cd + 201 
611 1.40 10.3.3 0.35 V.V.W. Cd + 2G) 
111 4.44 110 2.22 m. 
311 1.38 100 1.38 s. 
310 2.64 521 0.66 m. 
321 2.02 530 0.50 m.w. 
511 1.52 411 0.76 V.V.W. 
531 1.20 210 1.20 Ww. 
753 0.80 320 0.80 m.w. 
*v.s. = very strong, m = medium, w = weak, etc. 
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Structure factors calculated for structure I with u = 0.25 are included 
in table 1 for forms giving first-order reflections. (These forms are ar- 
ranged in order of decreasing interplanar distance.) These are in satis- 
factory agreement with intensities estimated from the reproduction of 
the photograph. Similar agreement is obtained with the powder photo- 
graph intensities for MgCh, ZnCk, CoCl, and NiCl, reported by Bruni 
and Ferrari.® 

The structure found for this series of divalent chlorides is shown in the 
figure. It is a layer structure closely related to that of cadmium iodide. 
Each cation is surrounded by six chloride ions approximately at the 
corners of a regular octahedron, six edges of which are shared with other 
octahedra in such a way as to form a layer. In the cadmium iodide 
structure these layers are superimposed in such a way that each cation 
is directly above a cation in the layer below, whereas in the cadmium 
chloride structure each cation is directly above an anion in the layer 
below, and below an anion in the layer above. It is probable that the 
cadmium chloride structure is stable for substances MX, in which the 
cation M has a coérdination number 6 and the anion has an appreciable 
but not great polarizability. 
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I propose to carry on the experimental investigation of some of these 
crystals in order to obtain accurate values for the parameter. The dis- 
cussion of interionic distances can be postponed until the results of this 
study are reported in the Zeitschrift fir Kristallographie. 





FIGURE I 


The structure of cadmium chloride and other bivalent chlorides. Large open circles 
represent chloride ions, small closed circles cadmium ions. Heavy lines mark the 
rhombohedral unit of structure, and light lines the coérdinated octahedron formed by 
six chloride ions about each cadmium ion. 


Summary.—With the aid of data obtained from a reproduction of a 
Laue photograph published by Bruni and Ferrari, it is shown that the 
unit of structure of CdCl. is a rhombohedron with a = 36° 40’ anda = 
6.35 A, about, containing 1CdCl, Cd at 000, 2Cl at uuu, fii, with 
u = 1/,. The isomorphous crystals MgCh, CoCh, and NiCk have a = 
33° 36’ and a = 6.21, 6.14 and 6.13 A, respectively. Other representa- 
tives of the structure are MnCl, ZnCl, RuCl, RhCl, PdCl, IrCle and 
PtCle. 

1G. Bruni and A. Ferrari, Rendiconti Acc. dei Lincei, 2, 457 (1925). 

2G. Bruni and A. Ferrari, Jbid., 4, 10 (1926). 

3 A. Ferrari, [bid., 6, 56 (1927). 

4V. M. Goldschmidt, Geochemische Verteilungsgesetze der Elemente, VIII, p. 147. 

5L. Pauling, J. Am. Chem. Soc., 49, 765 (1927). 

6 One line of small intensity, assigned indices {422} by Ferrari, is not compatible 
with the unit containing 1MClo. 
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THE ELASTIC PROPERTIES OF CERTAIN BASIC ROCKS AND 
OF THEIR CONSTITUENT MINERALS 


By L. H. ADAMS AND R. E. Gipson 
GEOPHYSICAL LABORATORY, CARNEGIE INSTITUTION OF WASHINGTON* 


Communicated August 2, 1929 


The rapid growth of the science of seismology, providing as it does more 
and more reliable data on the velocity of elastic waves at various depths 
within the Earth, calls for a closer study of the elastic constants of rocks 
and minerals because, given an accurate knowledge of the velocity of 
earthquake propagation at any depth and the elastic constants of all 
rocks under conditions prevailing at that depth, it is possible to speculate 
intelligently upon the nature of the rocks within the interior of the Earth. 

Several years ago it was shown! that measurements of compressibility 
taken with the value of Poisson’s ratio, which is approximately the same 
for all rocks, provided one of the best means of determining these elastic 
constants at high pressures. At that time the acidic rocks were investi- 
gated more thoroughly than the basic ones, and yielded more consistent re- 
sults. The importance which is now being attached to the properties of basic 
rocks in inquiries concerning the nature of the deeper parts of the Earth’s 
crust has made it seem advisable to carry out additional measurements of 
the compressibility of these rocks and certain of their constituent minerals. 

In this work we have measured the compressibility of several diabasest{ 
with an improved apparatus, and where necessary redetermined the 
compressibility of the minerals in order to formulate more clearly the rela- 
tion between the compressibility of the rocks and of their minerals. Fur- 
thermore, we have examined the compressibility of the garnets, grossu- 
larite and almandite, and of the pyroxene, jadeite, with a view to estimating 
the compressibility of eclogites, or rocks resembling basalt in gross com- 
position but crystallized as high pressure minerals. 

Experimental Method.—The apparatus, although entirely new, is very 
similar in construction to that already described in former publications! 
from this Laboratory. A leak-proof piston? is forced into a heavy-walled 
steel cylinder, which contains the specimen completely immersed in a 
suitable liquid (n-butyl ether). The hydrostatic pressure so generated is 
measured to 1 megabarye (about 1 atmosphere) by an electrical resistance 
gauge, and the corresponding decrease in volume is determined by measur- 
ing the travel of the piston with a dial-micrometer on which one division 
represents 0.01 mm. The apparatus is calibrated frequently by making 
a series of measurements on soft steel, the compressibility of which is 
accurately known.’ The formula for calculating the decrease in volume 
for each pressure has been given fully in previous papers.‘ 











714 GEOPHYSICS: ADAMS AND GIBSON Proc. N. A. S. 


Materials.—The sample of Whin Sill diabase was collected by Professor 
A. Holmes at Scordale Beck, near the Hilton lead mines, Westmorland, 
England. The rock is a fine-grained quartz-dolerite of density 2.937 at 
31.4°. It has been described by Holmes and Harwood,’ who give on 
page 512 of their paper a complete analysis. The calculated mode is ap- 
proximately as follows: quartz, 7; plagioclase (Ab,An), 47; augite, 19; 
hypersthene, 18; magnetite and ilmenite, 9. 


TABLE 1 
RESULTS FOR MINERALS 
ALMANDITE (243) ALMANDITE (241) GROSSULARITE (239) GROSSULARITE (245) 
01044 = —0.4 1044 = 0. 1044 = 1.2 1044 = 2.1 

—_— re 0.576 ri = 0.557 > = 0.600 Pd = 0.595 
(uzca- ~ V, * 1 ce, ~ ¥e *™ we. ~ Ve * 10 ave. 7. 1° oo. 
BARYES) (OBs.) (CALC.) X 104 (oBs.) (CALC.) X 104 (oBs.) (cALC.) X 104 (oBs.) (caLte.) xX 104 
12000 0.565 0.572 —0.7 0.550 0.561 —1.1 06.601 0.613 —1.2 0.606 0.616 —-1.0 
11000 0.517 0.515 0.2 0.506 0.505 0.1 0.550 0.553 -—0.3 0.556 0.556 0.0 
10000 0.455 0.457 —0.2 0.456 0.450 0.6 0.489 0.493 —0.4 0.494 0.497 —-0.3 
9090 0.404 0.399 0.5 0.397 0.394 0.3 0.435 0.433 0.2 0.437 0.437 0.0 
8000 0.347 0.342 0.5 0.349 0.338 1.1 0.382 0.373 0.9 0.372 0.378 —0.6 
7000 0.287 0.284 0.3 0.288 0.288 0.5 0.315 0.313 0.2 0.332 0.318 1.4 
6000 0.231 0.227 0.4 0.228 0.227 0.1 0.263 0.253 1.0 0.263 0.259 0.4 
5000 0.166 0.169 —0.3 0.165 0.171 —0.6 0.207 0.193 1.4 0.208 0.199 0.9 
4000 0.098 0.111 —1.3 0.117 0.116 0.1 0.141 0.132 0.9 0.161 0.140 2.1 
3000 0.055 0.054 0.1 0.054 0.060 —0.6 0.055 0.072 —1.7 (0.109)* (0.080)* ... 
2000 0.000 —0.004 0.4 0.000 0.004 -0.4 0.000 0.012 —1.2 0.000 0.021 —2.1 
GROSSULARITE (261)* JADEITE (251) JADEITE (249) LABRADORITE (235) LABRADORITE (237) 

10’a = 3.8 10‘a = 1.6 10’a = 3.1 10‘'a = 0.2 10‘a = 0.1 

10% = 0.611 10% = 0.735 10% = 0.775 10% = 1.457 10% = 1.417 

pa ont ps dome = 0.115 10K = —0.086 
7 i i ~ Y 1O coe. ~ Ve * ™ weve. ~ Ve * ™ wees, ~ Ve. * ™ wore. 


(ops.) (caLc.) X 104 (ons.) (caLe.) X 104 (ops.) (caLc.) X 104 (oBs.) (cauc.) X 104 (oBs.) (cALC.) X 104 


0.642 0.649 —0.7 0.748 0.751 —0.3 0.784 0.806 —2.2 1.342 1.344 —0.2 1.328 1.333 —0.5 
0.591 0.588 0.3 0.669 0.678 —0.9 0.715 0.729 —1.4 1.225 1.220 0.5 1.211 1.208 0.3 
0.524 0.527 —0.3 0.603 0.604 —0.1 0.659 0.651 0.8 1.101 1.094 0.7 1.092 1.081 1.1 
0.470 0.466 0.4 0.531 0.531 0.0 0.579 0.574 0.5 0.959 0.966 —0.7 0.946 0.952 —0.6 
0.409 0.405 0.4 0.460 0.457 0.3 0.508 0.496 1.2 0.830 0.835 —0.5 0.818 0.821 —0.3 
0.351 0.344 0.7 0.391 0.384 0.7 0.439 0.419 2.0 0.693 0.702 —0.9 0.683 0.689 —0.6 
0.276 0.282 —0.6 0.312 0.310 0.2 0.352 0.341 1.1 0.569 0.566 0.3 0.556 0.555 0.1 
0.217 0.221 —0.4 0.243 0.237 0.6 0.276 0.264 1.2 0.444 0.429 1.5 0.426 0.419 0.7 

0.178 0.163 1.5 0.192 0.186 0.6 0.286 0.289 —0.3 0.278 0.282 —0.4 

0.083 0.090 —0.7 0.104 0.109 —0.5 0.144 0.146 —0.2 0.146 0.143 0:3 

0.000 0.016 —1.6 0.000 0.031 —3.1 0.000 0.002 —0.2 0.000 0.001 —0.1 


* Omitted from the least square calculation. 


The Maryland diabase was collected by Dr. E. S. Shepherd on the 
Baltimore Pike near Frederick, Md. It is a fine-grained diabase, density 
3.033 at 31.0°. An unpublished analysis, by E. S. Shepherd, of the 
specimen used in this work is as follows: SiOQ2, 51.28; AipO2, 15.07; FesOs, 
1.12; FeO, 9.31; MgO, 7.97; CaO, 11.42; Na,O, 2.03; KO, 0.27; H:O +, 
0.39; H,O-, 0.09; TiOe, 0.78; P2O5, 0.13; Cl, 0.00; S, 0.06; Cr2O3, 0.05; 
MnO, 0.16. Calculated mode; plagioclase (Ab;Ane), 48; augite, 24; hyper- 
sthene, 25; magnetite and ilmenite, 3. 

The specimen of Sudbury diabase is the same as that used by Adams 
and Williamson.’ Density 3.002 at 31.1°. It is a moderately coarse- 
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grained olivine diabase. Calculated mode: plagioclase (Ab;Ane), 60; 
augite, 15; olivine, 17; magnetite and ilmenite, 8. 

The jadeite was from a large piece of massive material from Burma and 
was kindly furnished by Dr. W. F. Foshag of the U. S. National Museum. 
Density 3.328 at 30.9°. Microscopic analysis showed it to be practically 
pure jadeite. 

The fine specimen of translucent massive grossularite was also furnished 
by Dr. Foshag. Density 3.544 at 30.9°. An analysis of material from 
the locality has been made by Clarke and Steiger.’ 


TABLE 2 
RESULTS FOR DIBASE ROCKS 


MARYLAND DIABASE (227) 
104a 


MARYLAND DIABASE (229) 
1 


MARYLAND DIABASE (233) 





Ota = —0.2 4a = 10‘a = 
10% = 1.280 10% = 1.239 10% = 1.228 
= —0.115 ty = —0.135 ‘av = —0.130 
PRESSURE é 
(MEGA- oe 10 DIFF. °“ 1 DIFF. - IFF. 
BARYES) (oBs.) (caLc.) XX 104 (oBs.) (CALC.) x 104 eaik” (CALC.) xX 104 
12000 1.164 1.162 0.2 1.111 1.107 0.4 1.103 1.100 0.3 
11000 1.048 1.056 —0.8 1.007 1.009 —-0.2 1.005 1.002 0.3 
10000 0.958 0.948 1.0 0.911 0.908 0.3 0.903 0.901 0.2 
9000 0.835 0.837 —0.2 0.797 0.804 —-0.7 0.787 0.798 —1.1 
8000 0.722 0.724 —0.2 0.694 0.698 —0.4 0.683 0.692 —0.9 
7000 0.604 0.609 —0.5 0.582 0.589 -0.7 0.585 0.584 0.1 
6000 0.493 0.491 0.2 0.483 0.477 0.6 0.477 0.473 0.4 
5000 0.374 0.371 0.3 0.371 0.363 0.8 0.375 0.359 1.4 
4000 0.252 0.249 0.3 0.254 0.246 0.8 0.240 0.243 -0.3 
3000 0.118 0.125 —0.7 0.120 0.126 —0.6 0.119 0.124 —-0.5 
2000 0.000 —0.002 0.2 0.000 0.004 —-0.4 0.000 0.002 -—0.2 
SUDBURY DIABASE (217) SUDBURY DIAbasy (219) WHIN SILL DIABASE (221) WHIN SILL DIABASE (223). 
10’a = 0.1 10‘a 0 10‘a = 0.15 10‘a = 0.5 
10% = 1.440 10% = 1.351 10% = 1.712 10% = 1.688 
—> = — 0.130 10% = — 0.061 — = — 0.284 bay = — 0.280 
— Fe XO crs. Sethe , at ietwe a  : o 
(oBs.) (catc.) K 104 (oBs.) (catc.) X 104 (oBs.) (cALc.) X 104 (oBs.) (CaALC.) X 104 
1.317 1.312 0.5 1.295 1.291 0.4 1.434 1.430 0.4 1.418 1.414 0.4 
1.182 1.192 -1.0 1.163 1.168 -0.5 1.309 1.312 -0.3 1.299 1.298 0.1 
1.080 1.070 1.0 1.048 1.043 0.5 1.192 1.189 0.3 £188; 2 477 0.4 
0.941 0.946 -0.5 0.911 0.917 -0.6 1.054 1.061 -0.7 1.040 1.050 -—1.0 
0.816 0.819 -—-0.3 0.784 0.7909 -—-0.6 0.929 0.927 0.2 0.906 0.918 -—1.2 
0.683 0.689 -0.6 0.659 0.661 -—0.2 0.773 0.787 -—1.4 0.774 0.780 —0.6 
0.563 0.557 0.6 0.537 0.531 0.6 0.648 0.638 1.0 0.644 0.636 0.8 
0.428 0.422 0.6 0.407 0.400 0.7 0.498 0.490 0.8 0.502 0.487 1.5 
0.282 0.284 -—0.2 0.270 0.268 0.2 0.341 0.333 0.8 0.339 0.332 0.7 
0.143 0.144 -0.1 0.130 0.1385 -0.5 0.161 0.170 -0.9 0.166 0.172 -0.6 
0.000 0.001 —-0.1 0.000 0.001 -0.1 0.000 0.001 0.1 0.000 0.006 —0.6 


The almandite, from Salida, Col., was also furnished by Dr. Foshag. 
Density 4.160 at 30.9°. An analysis of almandite from the same locality 
has been made by Penfield and Sperry.® 

The specimen of labradorite, from Nain, Labrador, is the same as that 
described in a previous paper.’ The density is 2.696 at 29.5°, and the 
composition is approximately AbysAnge. 

Experimental Results.—The results for the three rocks and for the four 
minerals are shown in tables 1 and 2, which give the fractional decrease 
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in volume, — AlV’/V,,, caused by an increase of pressure from 2000 mega- 
baryes to the pressure indicated. Within the experimental error, the 
observations for jadeite and for the two garnets, fall on a straight line, 
when AV/V, is plotted against the pressure, P. Accordingly, they 
have been represented by the equation: 


— AV/V, = a+b (P — 2000) (1) 


where a and b are constants. These constants were calculated from the 
data by the method of least squares, and are shown in table 1 at the head 
of each series of results. The agreement between the observed values 
of AV /V, and those calculated by equation (1) is quite satisfactory and 
indicates that the readings of the dial-micrometer were consistent to 
better than 0.01 mm. 

For labradorite the graph of Al’/\’, against P is distinctly curved; 
the results, therefore, are represented by a quadratic equation: 


— AV/V, = a+b (P — 2000) + c (P — 2000)? (2) 


the constants, a, b, and c being obtained from a least-square solution. 

The compressibility, B(equal to —d\’//!’, dP) may be obtained by 
differentiating equation (1) or (2). The compressibility of labradorite 
varies considerably with pressure, and is found to be 1.48 x 10-* at 
P = 0, and 1.28 X 10~* at P = 10,000 megabaryes. These values are 
somewhat lower than those previously obtained (1.55 and 1.34, respec- 
tively), but the difference is partly accounted for by the lower value now 
taken for the compressibility of the reference substance, soft steel. 

The compressibilities of the minerals are summarized in table 3, in which 
have also been inserted the results previously obtained for other minerals 
commonly found in basic rocks and the observations on magnetite and 
ilmenite made by Madelung and Fuchs.® For those substances, whose 
change of compressibility has not been measured, the compressibilities 
at the various pressures have been determined by the empirical method 
already described.!° This method takes advantage of the fact that, 
approximately, the change of compressibility with pressure is a function 
of the compressibility alone. If the compressibility at only one pressure, 
or the mean compressibility over a certain range of pressure, is known, 
the compressibility at various pressures can be estimated, probably with 
considerable accuracy, and certainly with less error than if no allowance 
were made for the change of compressibility. 

In the plagioclase series only two members have been investigated, 
but on the assumption that the compressibility bears a linear relation 
to the volume-percentage of albite or anorthite the compressibilities of 
two other members of the series, AbgsAn4 and Abs;3Ang7, have been calcu- 
lated. In table 3 those values which have been obtained by interpolation, 
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extrapolation, or correction for the ‘‘statistical’”’ change of compressibility 
are enclosed in parentheses. For pressures from 0 to 15,000 megabaryes, 
this table gives a comprehensive view of the compressibilities of the more 
important minerals found in basic rocks. 

TABLE 3 


CoOMPRESSIBILITIES OF THE IMPORTANT MINERALS IN Basic Rocks 
COEFFS. OF EQUATIONS COMPRESSIBILITY, 1068, AT THE FOLLOWING 


(1) AND (2) PRESSURES (MEGABARYES) 
MINERAL 10% —10% 2000 7000 10, 15,000 
Quartz 2.61 0.216 2.70 2.61 2.40 2.27 ie 
Oligoclase Ab7gAnz: 1.70 0.105 .74 1.70 1.60 1.54 (1.44) 


1 
Andesine AbgoAnyo oe siete (1.59) (1.55) (1.45) (1.39) (1.29) 
Labradorite Ab:sAnse* 1.44 0.100 1.48 1.44 1.34 1.28 (1.18) 
Labradorite Ab33Ang7 ace Late (1.382) (1.30) (1.24) (1.20) (1.14) 


Actinolite 1.29 (1.87). (1.35) 1:20 . (8,25) 4.19) 
Diopside 1.06 (¥.13)- Gti) — 1.06 . (4:08). (6-98) 
Augite 1.01 (1.07) (1.05) 1.01 (0.98) (0.94) 
Enstatite 1.00 (1.06) (1.04) 1.00 (0.97) (0.93) 
Hypersthene 0.98 (1.04) (1.02) 0.98 (0.95) (0.91) 
Olivine 0.81 (0.85) (0.84) 0.81 (0.79) (0.76) 
Jadeite* 0.75 (0:78): 277) . 0:78 » 0.78)... O2D 
Garnet (Grossularite) * 0.60 (0.63) (0.62) 0.60 (0.59) (0.57) 
Garnet (Almandite) * Saar (0.60) (0.59) 0.57 (0.56) (0.54) 
Ilmenite Se er 0.56 (0.55) (0.53) (0.52) (0.50) 
Magnetite re Oe 0.55 (0.54) (0.52) (0.51) (0.49) 


* New determinations described in this paper. 


With each of the three rocks investigated the volume-decrease was not 
a linear function of the pressure, a fact which was at once evident from 
the graphs. Two series of measurements were made with the Whin Sill 
diabase, two with the Sudbury, and three with the Maryland. There 
seemed to be a tendency for the second or third series with a given speci- 
men to yield a somewhat lower compressibility than the first series. This 
was probably due to the liquid more completely filling the fine pores of the 
rock by the continued application of pressure. In figure 1 is shown graphi- 
cally the volume-change under pressure for the Whin Sill and for the 
Maryland diabase (two series for each). 

Discussion.—In a previous investigation’ the compressibility of rocks 
was measured under two different conditions, (1) the rock was placed 
directly in contact with the pressure-transmitting liquid, and (2) the rock 
was completely enclosed in a thin jacket of tin. The first method, leading 
to the ‘‘inner’’ compressibility, gave the volume changes due to the com- 
pression of the actual material particles of the rock because the pressure 
liquid had ample opportunity to fill up the interstices; but in the second 
method, from which the ‘‘outer’’ compressibility was calculated, the 
observed volume change included both the change due to the compression 
of the particles and that due to the closing of the pores. 
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If we are interested in the propagation of elastic waves through a rock 
of a certain composition and structure on the surface of the Earth the 
“outer’’ compressibility is the modulus which is of most importance to 
us, but in this paper we are more concerned with rocks formed under 
pressures at which a very compact structure would be most probable. 
For this reason and also on account of the greater ease of experimentation 
we have made all our measurements on unenclosed specimens. At lower 
pressures the differences between these two types of compressibility are 
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Final results for the volume change under pressure of Whin Sill and Maryland diabases. 


considerable for many, especially the more porous rocks, but beyond 
10,000 megabaryes the differences approach zero asymptotically and so 
our problem is simplified. 

In the paper by Adams and Williamson’ it is shown that at high pressures 
the compressibility of rocks tends toward the value obtained by averaging 
(according to volume-percentage) the compressibilities of the minerals 
present in the rock. At lower pressures, however, compressibility of the 
rock may be much greater than that calculated from the minerals and 
depends as much on the physical condition of the particular rock as on its 
mineral composition. Figure 2, which is taken (with minor changes) from 
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the earlier paper, illustrates the way in which the compressibilities of 


granitic and gabbroic rocks vary with pressure or depth. The large varia- 
tion at low pressures depends on the more or less accidental structure of the 
rock, while at high pressures the variation is mainly that due to differences 
in mineral content of the rock of the type in question. 

In table 4 are summarized, first the compositions of the various rocks 
in terms of the minerals present, and second the measured and calculated 
compressibilities at certain pressures. The results for the previously 
investigated Palisade diabase, recalculated to show the change of com- 
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Compressibility, @ x10® —> 


Gabbros and Diabases 





0 2000 4000 6000 8000 10000. ~=—-: 12000 
Pressure in Megaboryes —> 
FIGURE 2 


The compressibility of granites and diabases as a func- 
tion of pressure. The width of the shaded areas indicates 
the variation which may be expected for ordinary rocks 
within a given class. At high pressures the variation is 
governed by the mineral composition of the rock, at low 
pressures the variation in compactness is the dominating 
factor. 


pressibility with pressure, have been included in the table. There is 
also included the presumable mineral content and the calculated com- 
pressibility of the average plateau basalt as computed by R. A. Daly" 
from 50 analyses of Deccan, Oregonian and Thulean flows. 

The measured compressibility, which at 2000 mb. exhibits a consider- 
able variation from one rock to another, shows much less variation at 
10,000 mb. Moreover, at 10,000 mb. the compressibilities have approached 
fairly closely to those calculated from the minerals, being on the average 
only 0.10 X 10-* higher. In one case—the Maryland diabase—the figures 
are the same for the measured and calculated compressibilities. The 
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greatest divergence is found for the Palisade diabase, the measured and 
calculated values being 1.30 X 10~* and 1.16 X 10~, respectively. It is 
significant in this connection that by one criterion of the degree of weather- 
ing the Maryland diabase is a much fresher rock than the Palisade dia- 
base. According to unpublished measurements by E. S. Shepherd the 
latter yields 90 cm.* of gas (mostly water vapor) per gram when heated 
to a high temperature in a vacuum, while the former evolves only 30 cm.* 
per gram, which is close to the minimum observed for ordinary rocks. 
Only lavas from recent flows show less. In the former paper the presence 
of saussuritized feldspar in the Palisade diabase was noted as an indication 
TABLE 4 


COMPRESSIBILITY AND MINERAL COMPOSITION OF VARIOUS DIABASE ROCKS 
AVERAGE PLATEAU 


WHIN SILL SUDBURY PALISADE MARYLAND BASALT (DALY) 
MINERAL COMPOSITION 
Quartz 7 Sali ak mre? nea 
Plagioclase 47 (1,1) 60 (3,2) 54 (3,2) 48 (1,2) 50 (8,2) 
Augite 19 15 20 24 \ 37 
Hypersthene 18 Pa 14 25 j 
Olivine State 17 5 ws 4 
Magnetite and Ilmenite 9 8 7 3 9 
ComPRESSIBILITY, 8 X 10° 
Measured 1.70 1.37 1.54 1.23 
P = 2000 mb. 
Measured 1.26* 1.25* 1.30* OF * 
P = 10,000 mb. 
Calcd. from minerals 1.17 1.16 1.16 1.07 1.13 
P = 10,000 mb. 
Calcd. from minerals 1.09 1.08 1.09 1.02 1.06 
P = 15,000 mb. 


* These figures were calculated by passing a least-square straight line through the 
points from 8000 to 12,000 megabaryes. It is noteworthy that the same result was 
obtained from the slope of the quadratic equation at 10,000 megabaryes. 


that the rock was not perfectly fresh. The Maryland rock, on the other 
hand, is probably a very perfect specimen of a sound and fresh diabase. 
Even at 2000 megabaryes the measured compressibility is nearly as low 
as that calculated from the minerals (1.23 against 1.16). We may con- 
clude from the table that at high pressures the compressibilities of the 
various rocks approach and finally become equal to those calculated from 
the mineral content. Indeed, it appears that for any rock a minimum 
value of its compressibility at any pressure may be obtained by averaging 
the compressibilities of its constituent minerals. The actual compressi- 
bility may be considerably higher, especially at low pressures. It is in- 
teresting that granites lose their abnormality almost entirely at 2000 
mb. or less, but that the basic rocks in most instances still retain some of 
their abnormality even at 10,000 mb. 
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In a former paper a slightly higher value was assigned to the compressi- 
bility of diabase than will be found here. This is due to two factors: 
(1) better values now available for the compressibility of the standard of 
reference, steel, which are lower than those used in the previous work, 
and which exhibit an appreciable decrease with increasing pressure; and 
(2) the difference between the old and the new values for the compressi- 
bility of labradorite. 

From figure 2 it may be seen that the characteristic course of the com- 
pressibility of rocks as the pressure’is increased is first a rapid and then 
a much slower decrease of the compressibility. This course should properly 
be represented by an equation with an exponential term. Even though 
the volume change of a rock in the pressure-range from 2000 to 12,000 
megabaryes may, within the error of experiment, follow a parabolic equa- 
tion such as eqn. (2), it is very evident that an extrapolation by such an 
equation to either higher or lower pressures would be very unsafe. In- 
deed, an extrapolation in either direction would yield compressibilities 
lower than the true ones. At very low pressures the compressibility of a 
given sample of rock can be determined only by direct measurement on 
the particular sample, and at pressures above the experimental range a 
calculation based on the extrapolated compressibilities of the constituent 
minerals will give the most satisfactory value, although it is to be remem- 
bered that the value so obtained is merely the minimum compressibility. 
The bottom line of table 4 shows the calculated compressibilities at 15,000 
megabaryes of the four diabases and also of Daly’s average basalt (if 
crystalline). It does not seem possible that the compressibilities at the 
pressure stated can be significantly lower than the values given, and it 
does not seem probable that they can be very much higher. 

In a recent paper on ‘“The Outer Shells of the Earth,’”’ R. A. Daly!! has 
made a painstaking criticism of all previous determinations of the com- 
pressibility of rocks, and has concluded that the values which have been 
given for the compressibilities are too high, mainly because of two factors: 
first, the change of compressibility with pressure; and second, the ten- 
dency of elastic hysteresis to increase the effective elastic moduli of rocks. 
The first effect has already been discussed in a paper by the authors, and 
as for the second, we have not been able to observe any hysteresis in 
the volume change of igneous rocks when subjected to hydrostatic pressure. 

From the compressibilities given in table 4 the velocity of longitudinal 
waves through the various diabase rocks has been calculated by the 
equation:*® 

Uy = 13.13/+/10Bp, (3) 


where U, is the velocity in kilometers per second and p is the density. 
The results are shown in the following table. 
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TABLE 5 


VELOCITY OF LONGITUDINAL WAVES IN VARIOUS DIABASE RocKS 
VELOCITY IN KM./SEC. AT VARIOUS PRESSURES 
10,000 10,000 


, 15,000 

ROCK 8 MEASURED 8 MEASURED B CALC. B CALC. 
Whin Sill diabase 5.9 6.8 7.0 7.3 
Sudbury diabase 6.5 6.7 7.0 7.2 
Palisade diabase 6.1 6.6 7.0 7.2 
Maryland diabase 6.8 7.2 7.2 7.4 
Av. plateau basalt he | 7.3 


The velocities at 15,000 megabaryés equivalent to a depth of slightly 
more than 50 kilometers are appreciably less than 7.9 km. per sec., which 
is the velocity of longitudinal waves at this depth as deduced from seis- 
mologic data. It should be noted that the velocities given in table 5 
apply only if the rock is crystalline. The compressibilities of glass and 
glassy rocks‘ will be discussed in a subsequent paper. 

For plateau basalt the velocity at 15,000 megabaryes is 7.3 km./sec. 
Although the gap between 7.3 and 7.9 is not large, it seems difficult to 
close it and thus to allow crystalline basalt—with glassy basalt the gap 
is wider—at depths below the 50-60 km. discontinuity. We believe that 
the compressibility of diabase as here given cannot be in error by as 
much as 4 per cent, which would correspond to 2 per cent for the wave 
velocity. In order to have material of the composition of diabase below 
the 50-60 km. discontinuity it seems necessary to assume: (1) that that 
material is crystalline; (2) that the effect of temperature is negligible; 
and (3) that the seismologic value for the velocity at this depth is in error 
by at least 0.6 km./sec. 

The Compressibility of Eclogites——The compressibilities of grossularite 
and almandite, the only two garnets we could obtain as large cylinders, 
are given in table 1, and may be summarized by the statement that the 
average compressibility between 2000 and 12,000 megabaryes is 0.60 x 
10~* for grossularite and 0.57 x 10~* for almandite. Although the com- 
positions of these two garnets differ widely, their compressibilities are 
nearly the same and we are emboldened to place the average compressi- 
bility of pyrope between that of grossularite and almandite. 

We also determined the compressibility of jadeite. The results are 
in table 1 and indicate that the average compressibility of this mineral 
is 0.75 x 10~® reciprocal megabaryes. 

With the help of these results it is possible to calculate the approximate 
compressibility of an eclogite. In table 6 we have taken hypothetical 
eclogites whose mineral compositions are based on three of the Norwegian 
eclogites so carefully studied by Eskola.'*? The compressibility of each 
has been calculated at 7000 megabaryes from our data and thence the 
earthquake wave velocity through such rocks has been obtained in the 
usual way. It must be emphasized that although we have tried to keep 
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the compositions of the hypothetical rocks within the bounds of plausi- 
bility these calculations must be regarded as approximate and are intended 
merely to assign probable elastic constants to rocks composed of high- 
pressure minerals. The wave velocities in these eclogite rocks even at 
7000 megabaryes are about the same as those given below the 50-60 km. 
discontinuity by the seismic data, and hence our results do not conflict 
with the hypothesis" of an eclogite shell at that depth. Further discus- 
sion of this is postponed to some future occasion. ; 


TABLE 6 
CALCULATED COMPRESSIBILITIES AT 7000 MEGABARYES OF CERTAIN ECLOGITES AND 
THE VELOCITIES OF LONGITUDINAL WAVES THROUGH THEM 


PERCENTAGE BY VOLUME AVERAGE AVERAGE AVERAGE REFERENCE 
GARNET JADEITE DIOPSIDE p B Up 
50 30 20 bs Pe 0.73 8.0 1 
42 6 52 3.5 0.84 Oe’ s 2 
32 24 44 3.5 0.83 at 3 


1 Analogous to chloromelanite eclogite from Vanelvsdalen, Sdéndmgre, cf. Eskola, 
p. 31. 

2 Analogous to eclogite from Rgdhaugen, Almklovdalen, Séndm¢re, cf. Eskola, p. 27. 

3 Analogous to eclogite from Silden, Nordfjord, cf. Eskola, p. 37. 


Summary and Conclusions.—We have measured directly, at pressures 
between 2000 and 12,000 megabaryes, the cubic compressibility of the 
minerals, labradorite, jadeite, grossularite and almandite, and of diabase 
rocks from Sudbury, Canada; Frederick, Maryland; and the Whin Sill 
in the North of England. By combining these mineral results with those 
already known, we are able to give a table of the compressibilities, at 
various pressures, of all the important constituents of basic rocks, to 
compare the compressibilities of the rocks with those of their constituent 
minerals, and to conclude that the compressibility calculated from the 
mineral content gives a limit to which that of the rocks approaches at 
high pressure. Our results demonstrate, furthermore, that at 15,000 
megabaryes and 30°C. the maximum velocity of longitudinal elastic waves 
through rocks of basaltic composition and mode is 7.4 km. per sec. On 
the other hand, the compressibilities of garnet and jadeite, being found to 
be surprisingly low, lead to the conclusion that magmas of this composition, 
crystallized primarily as, or subsequently metamorphosed to eclogites, 
may transmit longitudinal waves with velocities exceeding 8 km. per sec. 
Bearing in mind that the pressures existing at 20, 40 and 60 km. below 
the surface of the earth are somewhat greater than 5, 10 and 15 thousand 
megabaryes, respectively, we may readily see the connection between 
these observations and deductions concerning the composition of the 
interior of the earth. Indeed, the existing data, in our opinion, place 
the possible components of the earth, below 60 km. and above the core, 
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in the following ascending order of probability: holocrystalline basalt, 
eclogite, peridotite. 

The authors are indebted to their colleagues, Drs. H. S. Washington, 
N. L. Bowen, J. W. Greig and E. S. Shepherd for many helpful sugges- 
tions, to Dr. W. F. Foshag for his aid in obtaining certain mineral speci- 
mens for this investigation, and to Professor A. Holmes for his kindness 
in collecting from the desired locality a sample of the Whin Sill diabase. 


* Papers from the Geophysical Laboratory, Carnegie Institution of Washington, 
No. 690. 

t Throughout this paper the word diabase is used in the American sense to indicate 
a holocrystalline basalt. It should be understood by British readers to mean dolerite. 

1 Adams, Williamson and Johnston, J. Am. Chem. Soc., 41, 1919 (12). Adams and 
Williamson, J. Franklin Inst., 195, 1923 (475). 

2P. W. Bridgman, Proc. Am. Acad. Arts & Sci., 49, 1914 (634). 

3 P. W. Bridgman, Proc. Am. Acad. Arts & Sct., 58, 1923 (166); see also Adams and 
Gibson, These PROCEEDINGS, 12, 1926 (275); Gerlands Beitr. z. Geophys., 15, 1926 (17). 

4 Adams and Gibson, Op. cit.; Adams and Williamson, Op. cit. 

5 Holmes and Harwood, Min. Mag., 21, 1928 (493). 

6 Adams and Williamson, Op. cit. 

7 Clarke and Steiger, U. S. Geol. Surv., Bull. No. 262, 1905 (72). 

8 Penfield and Sperry, Am. J. Sci., 32, 1886 (811). 

® Madelung and Fuchs, Ann. Physik, 65, 1921 (289-309). 

1. H. Adams, J. Wash. Acad. Sci., 17, 1927, 529. 

1R. A. Daly, Am. J. Sci., 15, 1928 (108). 

12 Pentti Eskola. Vidensk. Skrifter I. Mat. Nature. Klasse 1921, No. 8. 

13 For references to this hypothesis see Bowen, The Evolution of Igneous Rocks. 
Princeton Univ. Press, 1928. 


THE MINIMA OF.INDEFINITE QUATERNARY QUADRATIC 
FORMS 


By A. OPPENHEIM 
DEPARTMENT OF MATHEMATICS, UNIVERSITY OF CHICAGO 


Communicated June 28, 1929 


Let f = Soajnxjxp (ajz = ay;) be a quadratic form of rank n with real 
coefficients in ” variables x;, . . ., x,. The variables assume integer values 
only, the set 0, 0, ...excepted. One of the many problems associated 
‘with quadratic forms is that of finding the upper bound of L(f), the abso- 
lute lower bound of f, when f runs over the set of forms which have a 
given hessian a = \|a;||. We shall consider only-real indefinite quaternary 
forms. For indefinite binary and ternary forms, Markoff* has given two 
theorems, which are needed in the discussion of quaternary forms. By 
means of these theorems of Markoff, I have proved the following: 
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THEOREM. Let f be an indefinite quaternary quadratic form of hessian a. 
If the signature of f 1s —2 or +2 so that a 1s negative, then 


L‘(f) S —%a, (1) 
and, 1f equality holds, then f is equivalent to the form 
+ fy = L(x? — ay? —2 — P+ xy + x2 + X10). (2) 


If f is not equivalent to +f, then necessarily 
Lt < —3a. (3) 
If the signature of f is zero so that ais positive, then 
L4 S §a, (4) 
and, if equality holds, then f 1s equivalent to the form 
fi = L(x? — 9? — 2? + of + Qex + Qxy + xt + yt + 20). (5) 


If f is not equivalent to f,, then necessarily 


IA 
rit 


ES 5a. (6) 


In the present note I give an alternative proof I have discovered of the 
part of the theorem relating to forms of signature zero, but in this proof 
the inequality (6) is replaced by the cruder inequality 


16 9 y 
L® & 75a’. (7) 


A few preliminary remarks applying to forms in m variables will be 
necessary. ‘The adjoint of f is derived from f by the substitution 


fe) 
ay; = 52 = airy, 


where a is not zero so that the adjoint is 
= LAR 


where Aj, is the cofactor of a, in a. 

Let F be obtained from f by the transformation x; = }) >X;,Jj, of de- 
terminant +1, the /;, being integers. Then the transpose transformation 
Y; = oleiy, carries the adjoint G of F into the adjoint g of f. The quad- 
ratic form derived from F by putting s < m of the variables X equal to 
zero is called a section of f. If in particular },, = 1 and J; = 0 (¢ = 2 

. ., ”) so that the coefficient of x? is unaltered, and we do not put X; = 
0, the section so obtained may be called a section through the axis of x, 
or simply a principal section with respect to %;. 
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Suppose now that f is an indefinite non-nul’ quaternary quadratic 
form of signature zero. Every ternary section of f is necessarily indefinite. 
For suppose the section is obtained by putting X, = 0. We may write 


f=bet+.. Se bz, 


where 2; is a linear function of X;, X; +4 ,... Since the signature of f 
is zero, two of the }; are positive and two are negative. And it follows 
at once that the ternary obtained by putting \, or what is the same 
thing, 2; equal to zero, is indefinite. Since a form and its adjoint have 
the same signature, any ternary section of g is also indefinite. 

Consider the principal ternary section F(X,, X2, X3, 0) of f under x; 
= Solj,X;, where i; = 1, lx = 0 (i = 2, 3, 4), of determinant +1. Its 
hessian D is the coefficient of Yj in G. But the transformation Y; = 
Dleiye carries G into g. It follows that D is represented by the ternary 
g(0, ye, Vs, ys) for the integers ye, ¥3, ys determined from the consistent 
linear equations obtained by putting », = Y; = Yo = Y3 = O and Y, 
= ], 

Conversely, if D is any number properly represented by the ternary 
2(0, ve, . . .), we can find a section of f through the axis of x, with hessian D. 

Suppose now that L is not zero and is attained so that it is represented 
and properly represented by | f |. We may therefore take a, = +L. 
By Markoff’s theorem on indefinite ternary forms, we have 

lan |?=3%|D| orelse | an |* S$ ?| D|, (8) 


1 


for the absolute lower bound of F(X), X2, X3, 0) is plainly L. 
Let M be the lower bound of the indefinite ternary g(0, ye, 3, v4) of 
hessian a),a”. Then 


M? = 3|ana?| or M?* S =| ana? |, (9) 
and, if M is not attained, then (9.2) holds with < instead of < or =. 
Since D is any number properly represented by this ternary, we may 
subject D to the conditions 


" 9 1 | | 


Di = $|ay,a?| orelse | D |3 <? | aya” |. (10) 


From (8) and (10) we see that L' < #% a, unless both (8.1) and 


(10.1) hold, in which case L* = § a. Employing Markoff's result that 
an indefinite ternary form for which the cube of the absolute lower bound 
is equal to two-thirds of the hessian must be equivalent to the appropriate 
multiple of the form x? — y* — z* + xy + xz, which has hessian 3/2 and 
minimum unity, we find that the forms of lower bound L which satisfy 
(8.1) and (10.1) must be equivalent to the form (5), wherein L‘* = §a. 

If L is not attained, rather more argument is required. In this case, 
if e is any positive number, however small, we can choose a, so that 
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| | ay | < Lil + e). Let N be the absolute lower bound of F(X), 
Xe, X3, 0). It is easy to see that our previous arguments will yield (1) 
unless both 


N*=3|D| and |D|* = §| ana? | 


for all sufficiently small positive e. It is plain that now L < N and that 
| au | = N We deduce that N‘ = 4a | 9 and hence that 


Lt < 4a|9 < L41 + e)4. 


But clearly this is impossible if » is small enough, since L and a are in- 
dependent of e. Our theorem is therefore proved. 

This method fails in the case of forms of signature + 2, since the ternary 
sections employed may be definite or indefinite, and the information de- 
rived is insufficient to determine the precise upper bound. 

The principle of the method is due to Korkine and Zolotareff* who ap- 
plied it to obtain the precise upper bound for the minima of definite 
forms in four variables from the known upper bound for definite forms in 
three variables. 

We observe finally that for indefinite forms in five or more variable 
with commensurable coefficients, L(f) must be zero, or, what is here the 
same thing, f is a nul form. It is very likely true that L(f) must be zero 
when the coefficients are incommensurable. But this has not yet been 


proved. 

* A. Markoff.!:3\ For binary forms see also F. G. Frobenius,‘® I. Schur,®, 
and R. Remak.? Markoff proves his ternary theorem only for forms with commensur- 
able coefficients. It is not difficult however to extend his analysis so as to cover the 
case of incommensurable coefficients, including therein the case when the lower bound 
is not necessarily attained. The ternary theorem has been extended by L. E. Dickson.* 

t If f is nul, our theorem is trivial. 

1A. Markoff, “Sur les formes quadratiques binaires indefinies,’’ Math. Ann., 15, 
381-409, 1879. 

2 Ibid., 17, 379-399, 1880. 

3 “Sur les formes quadratiques ternaires indefinies,’’ Jbid., 56, 233-251, 1903. 

4G, Frobenius, “Uber die Markoffsche Zahlen,” Sits. Preuss. Akad. Wiss., 458-487, 
1913. 

5 “Uber die Reduction der indefiniten biniren quadratischen Formen,”’ Jbid., 202- 
211, 1913. 

6]. Schur, “Zur Theorie der indefiniten binéren quadratischen Formen,” Jbid., 
212-231, 1913. 

7R. Remak, ‘‘Uber indefiniten bindre quadratische Minimalformen,” Math. Ann., 
92, 155-182, 1924. 

81. E. Dickson, Studies in the Theory of Numbers, University of Chicago, 1929. 

9 A. Korkine and G. Zolotareff, ‘‘Sur les formes quadratiques positives quaternaires,” 
Math. Ann., 5, 581-583, 1872; Korkine, Coll. Papers, 1, St. Petersburg, 285-288, 1911 
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ON THE NUMBER OF CYCLIC SUBGROUPS OF A GROUP 
By G. A. MILLER 
DEPARTMENT OF MATHEMATICS, UNIVERSITY OF ILLINOIS 


Communicated July 25, 1929 


If d;, dz, . . ., d, represent all the different positive divisors of the order 
g of a finite group G, and if ¢(d;) and a; represent, respectively, the totient 
of d; and the number of the cyclic subgroups of order d;,7 = 1, 2,...,X, 
contained in G, then it results from the fact that two distinct subgroups 
of the same order cannot have in common any generator of these sub- 
groups and be generated thereby that 


a,(d,) + aep(de) +... + ayd(d,) = g. 


It will be assumed that d, = 1 and that d, = g. Hence it results that 
a, is always unity while a, is either unity or zero. A necessary and 
sufficient condition that a, = 1 is that G is cyclic, and whenever a,_ ; 
then ag; ~@=...=Q 1. As this case is so elementary it will be 
assumed in what follows that a, = 0, and it will be seen that whenever 
this condition is satisfied at least one of the coefficients a), do, ..., a — 1 
is as large as 3. 

When G contains more than one Sylow subgroup of a given order this 
theorem follows directly from the fact that each of these Sylow subgroups 
contains at least one cyclic subgroup which is not transformed into itself 
by a cyclic subgroup of the same order contained in an arbitrary one of 
the others, since a Sylow subgroup cannot be transformed into itself by 
another subgroup of the same order. Hence it may be assumed that G 
is either a prime power group or the direct product of such groups if the 
non-cyclic group G does not contain at least three cyclic subgroups of the 
same order. In both cases G must contain a non-cyclic prime power group, 
and hence it remains only to note that such a group contains at least three 
cyclic subgroups of some order. When this prime number ? is odd it is 
known that the number of the subgroups of order p contained in G is of 
the form 1 + kp, where k > 0. When » = 2 and G contains only one 
subgroup of order 2 then the number of its cyclic subgroups of order 4 
is known to be of the form 1 + 2k, where k > 0, since G must be dicyclic 
in this case. Hence we have proved the following theorem: Every non- 
cyclic group contains at least three cyclic subgroups of some order. 

1. Groups containing at least one cyclic subgroup whose order is an 
arbitrary proper divisor of the order of the group. 

If a non-cyclic group G of order g has the property that it contains at 
least one cyclic subgroup whose order is an arbitrary proper divisor of g, 
that is, any divisor of g except g itself, then its order cannot be divisible 
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by more than two distinct prime numbers. ‘To prove this fact we may 
first note that if g would involve more than two such factors it may be 
assumed that q is the smallest of them and that G contains a cyclic sub- 
group H of order g/g. It is easy to prove that H must be an invariant 
subgroup of G. This results from the facts that every Sylow subgroup 
of G except the Sylow subgroups of order g” must be contained in H, 
and that q”~* of the operators of each of the latter Sylow subgroups must 
also appear therein according to well-known theorems. Since an operator 
s of order g” must appear in cyclic groups whose separate orders are the 
product of g” and the order of each of the other Sylow subgroups it results 
that s must be commutative with every operator of H. This is impossible 
since G is non-cyclic and hence it has been proved that g cannot be divisible 
by more than two distinct prime numbers. 

It may now be assumed that g = p”q” and that the invariant sub- 
group H is of order p”q"~'. Since the group of isomorphisms of the 
cyclic group of order p” is cyclic and G is non-abelian it results that p — | 
is divisible by g. If m > 1, then G, therefore, cannot contain a cyclic 
subgroup of order p”~'g". Hence m = 1 while n can have any arbitrary 
value greater than 0. Moreover, the group G is completely determined 
by this value of ~ whenever p and gq are given prime numbers. Hence 
the following theorem has been established: Jf a non-cyclic group G of 
order g has the property that 1t involves at least one cyclic subgroup whose 
order 1s an arbitrary proper divisor of g, then g 1s divisible by at most two 
distinct prime numbers p and q, and when g is divisible by both of these num- 
bers and p > q then p — 1 ts divisible by q and g = pq". Moreover, there 
is one and only one such group for every positive value of n and it contains 
only cyclic proper subgroups. Every such group contains p such subgroups 
of order ” but only one such subgroup of each other order. 

In order to complete the determination of all the non-cyclic groups 
which have the property that each of them contains at least one cyclic 
subgroup whose order is an arbitrary proper divisor of the order of the 
group it remains only to consider the groups of order »”, p being a prime 
number, which have this property. All of these groups are known since 
all the groups of order p” which contain-a cyclic subgroup of order p”~' 
have been determined. Hence there is one and only one such abelian 
group for every value of m > 1, and there is also one and only one such 
non-abelian group which is conformal with this abelian group whenever 
m > 3. When p is odd the latter group exists also when m = 3 and 


there is no other group which satisfies these conditions. When p = 2 
and m = 3 the octic and the quaternion group satisfy these conditions, 
while when p = 2 and m > 3 there are five groups of order 2” which 


satisfy these conditions. By combining these known results with the 
theorem closing the preceding paragraph we therefore obtain all the non- 
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cyclic groups which satisfy the condition that each group contains at 
least one cyclic subgroup whose order is an arbitrary proper divisor of 
the order of the group. 

2. Groups containing at most three cyclic subgroups of the same order. 
It was noted above that every non-cyclic group contains at least one set 
of cyclic subgroups of the same order which includes as many as three 
such subgroups. We proceed to consider all those non-cyclic groups 
which contain no more than three cyclic subgroups of the same order. 
Since every non-cyclic group of order p”, p being any odd prime number, 
contains at least 1 + p subgroups of order p it results that every Sylow 
subgroup of odd prime power order which appears in such a group G of 
order g must be cyclic. If g = 2” and contains more than one cyclic 
subgroup of order 2, a > 2, it must contain an even number of such sub- 
groups and hence this number must be 2 in the present case. If G con- 
tains exactly three cyclic groups of order 4 it is known to be either the 
quaternion group or the group of order 16 obtained by extending the 
cyclic group of order 8 by operators which transform into their third 
power each operator of this cyclic group. As the latter contains five 
subgroups of order 2 it does not: satisfy the conditions under considera- 
tion and the quaternion group is therefore the only group of order 2” 
which satisfies these conditions and contains exactly three cyclic sub- 
groups of order 4. Hence the quaternion group is the only group of order 
2" which satisfies the two conditions that it contains no more than three 
cyclic subgroups of any order and that it contains exactly three cyclic sub- 
groups of order 4. 

In order to complete the determination of all the non-cyclic groups of 
order 2”, which contain at most three cyclic subgroups of the same order, 
it remains only to consider those which have three subgroups of order 2 
and two cyclic subgroups of each of the orders 4, 8, . . ., ge". ie fact, 
since a, = 0 and ¢(2”) = 2”~' each of the literal coefficients in the 
equation 


$(1) + 36(2) + as(ds) +... + a-16(d,-1) = 2” 


must be equal to 2. Hence there is one and only one abelian, and one and 
only one non-abelian group of order 2”, m > 3, which satisfies these 
conditions. These two infinite systems of conformal groups together 
with the non-cyclic group of order 4, the quaternion group, and the abelian 
group of order 8 and of type (2, 1) are therefore the only non-cyclic groups 
of order 2” which satisfy the condition that each of them contains no 
set of more than three cyclic subgroups of the same order. 

If g is divisible by more than one. prime number it may be the direct 
product of any one of the groups just determined and an arbitrary cyclic 
group of odd order or it may contain more than one Sylow subgroup of 
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order 2”. In the latter case each of its Sylow subgroups must be cyclic 
and G must involve an invariant cyclic subgroup whose order is the direct 
product of the cyclic group of order 3.2”~' and an arbitrary cyclic group 
whose order is prime to 6. Moreover, any such direct product can be 
extended by an operator of order 2”, which transforms the operators of 
order 3 in this cyclic group into their inverses but is commutative with 
each of the other operators contained therein so as to obtain a group which 
has three cyclic subgroups, when the order of such a subgroup is 2”, into 
any divisor of the order of the group which is prime to 6, but only one 
such subgroup when its order is not divisible by 2”. Hence the following 
theorem: Jf a group contains no more than three cyclic subgroups of the 
same order and is not a prime power group nor the direct product of such 
groups it can be constructed by extending the group which ts the direct product 
of the cyclic group of order 3.2"~' and any cyclic group whose order is 
prime to 6 by means of an operator of order 2” which has its square in this 
direct product and transforms into their inverses the operators of order 3 
contained therein but 1s commutative with all its other operators. 
! Miller, Blichfeldt, Dickson, Finite Groups, p. 131, 1916. 


ON NEW INTEGRAL ADDITION THEOREMS FOR BESSEL 
FUNCTIONS'! 
By HENRY P. THIELMAN 
Onto STATE UNIVERSITY, CoLUMBUS, OHIO 


Communicated August 5, 1929 


This paper contains some of the main results obtained by the author, 
which together with their derivations are to be published elsewhere. In 
obtaining these results use has been made of the theory of functional! 
composition as developed by Vito Volterra? and Joseph Pérés.* It must 
be said that Hadamard‘ has derived certain integral addition theorems 
for Bessel Functions. His results are, however, particular cases of a 
more general theorem stated here. 

THEOREM. Two Bessel functions of any orders a > —1, B > —1 have 
the additive properties expressed by the following equations: 


y a oe 8 Pi eit 
fate — pt7.eVae HI — »Ppl2VEE— 9) ld = 


+6+1 


a [(a + d(x — »)]| 2 








Ja+6 +1127 (a + b)(x — y)] (1) 


(a+ 5)" ter! 





SPC 











732 MATHEMATICS: H. P. THIELMAN Proc. N. A. S. 


y a 1 
fia — *)PI[2Va(x — t)] [b(t — y)] 2A[2V ot — y) |dt = 


“[(a + bx — »)P 
ee +i a Jal2/(a + Oe — 9] — 








— ; 
wi owe Tl2/ale—y] (2) 


fia -~° 5 3y,[2V ale — — t)|[b(t — 9) )eypl2V 0 — 0 — y)|dt 
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_ Via + )@ — »)P 


ae Ig[2V (a + b)(x — 9)] 





e+ — »)P 


a 


Jg(2V d(x — y)] (3) 


f ele = I-37, (2V ace — A] [b(t — ») | nV be — y) dt 
_ (a + b/)l2/G@+ He — I 


ab[(a + b)(x — )P 
_ Al2Vae =] _ A2Vbe = 9] 


blals — »)F alls — y)P 








(4) 





To obtain Hadamard’s result from this theorem we proceed as follows: 
We write 
ndis i —- 
ip(X) = ( — X)2J,[2V — X] (5) 
and note that 
d[j,(bt)] = bjg—1(bt)dt q> 0. (6) 


If in (1) we interchange a and b, and seta =q-—1,8=p,x =0,y = 1, 
this equation can be written, by means of (5) and (6), in the form 


ae sags a’b? 
$2: jpla(1 — t)|d[j,(bt)] = G@+bPr Iptg(a + 0), 


which is the general theorem of Hadamard. 
If in (3) we interchange a and b, and set 8 = p, x = 0, y = 1, then 
noting that 


d[j.(bt)| = 


dt, 





- jy(bt) 
b 
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we obtain the equation 


4 : ris(a + b) 
is jpla(1 — #)|d[y,(bt)| = eS 


which is the exceptional case in Hadamard’s general theorem. 

We thus see that Hadamard’s result. corresponds to particular cases of 
equation (1) and (3), while there is nothing in his work which might cor- 
respond to equation (4). 

By means of a new theory of composition on functions in 2 variables 
(n = 1, 2, 3, .. .) the author has derived integral addition theorems for 
hypergeometric series of a certain generalized type, which correspond 
to equations (1) to (4) of this paper. 

The writing of this paper was suggested by Professor A. D. Michal. 

2 Vito Volterra, Lecons sur les fonctions de lignes, Paris (1913). The Theoy of 
Permutable Functions, Princeton University Press (1915). 

3 Vito Volterra et Joseph Pérés, ‘“‘Lecons sur la composition et les fonctions per- 
mutable.” 

4 Hadamard, “Sur un probleme mixte aux derivées partielles,’’ Bulletin de la Soctéte 
Mathématique de France, tome 31 (1903), p. 220. 
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DWARF, A NEW MENDELIAN RECESSIVE CHARACTER 
OF THE HOUSE MOUSE 


By GeorcE D. SNELL 
Bussey INSTITUTION, HARVARD UNIVERSITY 


Communicated July 22, 1922 


During the course of the past winter a new recessive mendelian character, 
dwarf, was found among the mouse stocks of the writer at the Bussey 
Institution. This is the second case of hereditary dwarfism to appear in 
rodents. The first was found by Sollas* who has reported a strain of 
guinea-pigs in which '/, of the young (25 out of 99):from certain matings 
were undersize and malformed at birth. ~Most of these dwarf guinea-pigs 
died very early, and the few that lived were sterile and did not attain 
more than about one half normal size. The shape of the skull and various 
bones of the body was also altered. 

In the case of the dwarf mice, mature individuals are only about one 
fourth the weight of their normal brothers and sisters, scarcely bigger in 
fact than an ordinary mouse 16 or 17 days old. It is not until the 14th 
day, however, that their reduced size begins to manifest itself. Up to 
this time they grow as rapidly as their unaffected sibs, though they can 
usually be identified on the 12th or 13th day by their shorter noses and 
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tails. After 14 days they grow scarcely at all, though aside from per- 
manent shortness of their noses, most parts of the body especially ears 
and tail, develop enough to make the proportions of the animal as a whole 
those of an adult. The opening of the eyes is delayed by 2 or more days 
beyond its usual time of occurrence at 14 days. 

Despite their greatly reduced size, the dwarf mice are fairly healthy. 
All that have so far appeared (and not intentionally been killed) are still 
alive, the oldest being now 5 months of age. In vigor they are some- 
what sub-normal, showing very little tendency to run about and never 
trying to jump out of the jar in which they are placed during examina- 
tion. Both males and females, so far as can be judged from tests made 
to date, are entirely sterile. 

The character was derived from a strain of black silver mice very 
kindly given to the writer by Dr. L. C. Dunn, who obtained them from Wm. 
Turton, a fancier of Ilkeston, England. ‘They originally appeared in the 
F; of crosses of a silver male, S13, with two different pink-eyed females 
and also among the offspring of a silver male and female mated together. 
Some of the normal F2s from male S13 have now been proved to be hetero- 
zygous for the character, as a number of dwarf mice have appeared in 
F;, and in litters from a backcross of F; to F;. Matings of heterozygotes 
have so far yielded 103 normals to 34 dwarfs, a very close approximation 
to the expected 3:1 ratio. Owing to the sterility of the dwarfs, matings 
of homozygous recessives with heterozygotes and inter se have not been 
possible, but the clear-cut nature of the character, its appearance in a 
3:1 ratio in F:, and the fact that it has now been carried through two 
generations in the heterozygous condition since its first appearance, clearly 
establish it as a recessive unit factor. 

Linkage tests are now in progress and will be reported in detail later, 
but it is already clear that dwarf is not sex-linked and probably not linked 
with pink-eye or dilution. 

It is planned to undertake a physiological study of the dwarf mice as 
it seems likely that they will furnish excellent experimental material for 
various problems in connection with growth. 


* Sollas, I. B. J., ‘Inheritance of Colour and of Supernumerary Mamme in Guinea- 
Pigs, with a Note on the Occurrence of a Dwarf Form,” Reports Evol. Com. Royal Soc., 
5, 51-79 (1909). 
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CRITERIA OF MALE AND FEMALE IN BREAD MOULDS 
(MUCORS)' 


By SopHiA SATINA AND A. F. BLAKESLEE 


CARNEGIE INSTITUTION OF WASHINGTON, DEPARTMENT OF GENETICS, COLD SPRING 
Harsor, N. Y. 


Communicated August 10, 1929 


The Mucors, commonly called bread moulds, are a group of fungi 
apparently simple in structure and in their methods of reproduction. 
All are able to multiply rapidly, perhaps indefinitely, by means of non- 
sexual spores. We have still in cultivation, for example, the sexual 
races of the common Rhizopus which were isolated in 1904. A few 
species are hermaphroditic (homothallic) and produce sexual spores 
(zygospores) through the interaction of branches from the same individual. 
The great majority, however, are dioecious (heterothallic) and require 
for sexual reproduction the interaction of individuals of the two opposite 
sexes provisionally designated as (+) and (—). Which of these two 
sexes is male and which is female is the problem of our present discussion. 

The process of zygospore formation in dioecious species can be illustrated 
by figure 1, a-e. Two filaments of opposite sex (a) produce, at the point 


bi (+4 (—) (+) (-) 
(a) (0) (c) (d) (e) 


FIGURE 1 


of contact, swellings which enlarge (b), and from which cross walls 
cut off gametangia (c) which are morphologically indistinguishable. By 
dissolution of the intervening double wall, the contents of the gametangia 
are brought in contact (d) and the fusion cell thus formed develops into 
a mature thick-walled zygospore (e) which apparently functions as a 
resting spore. The cells of Mucors are multinucleate, and the filaments 
are usually without cross walls except at the formation of the reproductive 
bodies. Since the uniting sex cells contain more than a single nucleus 
they are called gametangia instead of gametes. 

The nuclei are minute and their behavior is difficult to follow. Prob- 
ably for this reason, although there has been considerable cytological 
work done on Mucors, there is no general agreement regarding the cyto- 
logical happenings during sexual reproduction. Few? have ventured to 
suggest even the presence of chromosomes. It is unreasonable, therefore, 
to expect much aid from a study of the chromosomes in determining which 
of the two sexes in Mucors corresponds to the male and which to the 
female of higher types of plants. 
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Tests with over 2000 races included in thirty-four or more species and 
twelve genera furnish ample evidence*® that in the dioecious Mucors we 
have two and only two sexes present comparable to those in higher forms. 
Races of given species were grown in pairs under conditions suitable to 
zygospore formation and the reactions recorded. Between 10,000 and 
20,000 such intraspecific combinations have been made without evidence 
for sex intergrades. The supposed examples of sex intergrades in the 
Mucors have been based upon relatively meager data and, as we have 
shown,‘ were probably wrongly interpreted. A given race shows itself 
to be either (+), (—) or fails to show any sexual reaction and is classed 
as a “neutral.’”’ Such neutral races are perhaps actually (+) or (—) 
but not of sufficient sexual activity to take part in zygospore formation 
with the races of the opposite sex 
with which they have been tested. 
Various grades of sexual vigor are 
apparent in the different races tested 
within a single species and some would 
have been classed as neutrals if they 
had not been tested with a particular 
race of the opposite sex with which 
they gave a scanty production of zygo- 
spores. We have found no evidence 
for ‘relative sexuality” in the dioecious 
Mucors such as Hartmann’ has de- 
scribed in Ectocarpus and_ believes 

FIGURE 2 exists in a wide range of forms. 

Not only are there only two sexes 
within a given species but the same two sexes are present throughout all 
the species tested as shown by the “imperfect hybridization’ reaction. 
Figure 2 is a diagram of a culture in which zygospores are represented 
by circles at the line of contact between (+) and (—) races of species V. 
In the same culture were planted the opposite sexes of another species IV. 
The dotted lines represent the imperfect sexual reactions which take place 
when opposite sexes belonging to two different species come in contact. 
Such imperfect reactions may be so abundant as to be visible as a white 
line. The reactions resemble early arrested stages in zygospore formation 
and never reach a more advanced stage than the formation of gametangia 
shown in figure Ic. 

Imperfect hybridization reactions have been employed in determining 
the sex of unmated races® and have led to the conclusion that there is 
something fundamental in common to the same sex throughout all the 
species in the Mucors. Imperfect reactions have also shown that the 
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same (+) and (—) sexes exist in the single individuals of hermaphroditic 
species.’ 

Although no constant morphological differences have been observed 
in the uniting gametangia of dioecious species nor in those of several 
hermaphroditic species, there are certain species of the latter in which 
there is a striking and constant difference in gametangia size, connected 
with differences in the hyphae which bear them. Since in practically all 
forms in which there is a size difference between gametes (the larger is 


Dioecious Species 
(+) (-) 














Hermaphroditic Species 


Zygorhynchus heterogamus Mbsidia spinosa 


FIGURE 3 


the female) it was thought possible by imperfect reactions to determine 
which sex of dioecious Mucors showed an imperfect reaction with the 
larger gametangia of these heterogamic hermaphrodites and thus to be 
able to say which of the (+) and (—) sexes was male and which female. 
Earlier investigations showed that the (+) race reacted with the smaller 
gametangium of three species of heterogamic hermaphrodites and the 
(—) race reacted with the larger gametangium of one of the three. These 
reactions would seem to warrant one in calling the (+) sex female and 
the (—) male. A preliminary test of a fourth heterogamic hermaphrodite, 
however, seemed to give somewhat different reactions with (+) and (—) 
testers of dioecious species and rendered it inadvisable to draw definite 











738 GENETICS: SATINA AND BLAKESLEE Proc. N. A. S. 


conclusions as yet as to the homologies of sexes in Mucors and higher 
forms. The matter seemed worthy of a more detailed study, the results 
of which are summarized in the succeeding paragraphs. 

Figure 3 presents diagrams of the sexual reactions between three types. 
The one above is a dioecious species used as a standard tester and its 
gametangia are labeled (+) and (—). The two types below, Absidia 
spinosa on the right, Zygorhynchus heterogamus on the left, are her- 
maphrodites with unequal gametangia, the larger in each case borne on 
the termination of an enlarged branch. 

The connecting lines in the diagram show the imperfect sexual reactions 
which have been observed. Thus (+) races of the dioecious species 
give reactions with the small gametangium of Absidia spinosa showing 
that the smaller gametangium of this species is (—) while the (—) race 
reacts with the larger gametangium showing that the latter is (+). A 
number of other heterogamic hermaphrodites (Dicranophora, Zygorhyn- 
chus Moelleri, Z. Vuillemini, etc.) give imperfect sexual reactions like 
those of Absidia spinosa. Zygorhynchus heterogamus, however, of which 
we have tested 4 races behaves differently from any other species investi- 
gated. In this species the smaller gametangium from its reactions must 
be considered (+) and the larger (—), a condition exactly the reverse 
of that found in Absidia spinosa. Imperfect sexual reactions occur 
between two hermaphrodites. When Absidia spinosa is grown with 
Zygorhynchus heterogamus, the large gametangium of the one species 
reacts with the large gametangium of the other while the small game- 
tangia of the two species react similarly inter se. The imperfect sexual 
reactions between these two hermaphrodites are such as might be expected 
from the nature of their gametangia as determined by their reactions with 
(+) and (—) races of dioecious species. 

If, in the evolution of these forms, the heterogamic condition has been 
derived from isogamy, it is obvious that in one species the (+) has de- 
veloped into the small gametangium while in others the (+) has developed 
into the large gametangium. It is as if in some species the eggs were 
male and the sperms female sex cells. Such a condition, so far as we 
are aware, is without parallel. It raises the question whether in groups, 
within which the sexual interactions have not been so thoroughly stid- 
ied, it is safe to homologize the sexes on the basis of size of conjugating 
gametes. In one of the flagellates, Mastigella vitrea, according to Wilson,® 
Goldschmidt found that the male gamete is non-motile while the female 
gamete is flagellated and in Monocystis sp., one of the Gregarines, accord- 
ing to Hartmann,® Naville found that the male are larger than the female 
gametes, but the differences shown in the figures given are at best slight. 

In no case in the literature, however, have we found record of the male 
cell being constantly so markedly superior in size to the female as is the 
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case with at least one of our heterogamic hermaphroditic Mucors. This 
is true whether the (+) is female and the (—) male as the reactions with 
Absidia spinosa would suggest or vice versa as the reactions with Zygor- 
hynchus heterogamus alone might indicate. It would seem that in the 
Mucors the terms (+) and (—) are labels of something more fundamental 
than the terms male and female if these latter are defined by size differences 
in the conjugating cells. 

In our discussion we have applied the sex terms (+) and (—) to game- 
tangia. ‘There has been a proper objection to applying the terms male 
and female to sperms and eggs in flies and mammals without regard to 
their chromosomal constitution, since half of the sperms have the same 
sex chromosomes as the eggs. It appears that the reactions of these 
gametes is not determined by their chromosomal constitution alone. 
Similarly the reactions of what has been convenient to call (+) and (—) 
gametangia in Mucors may not be due to differences in chromosomal 
content. They may differ in somewhat the same way as do sperms and 
eggs. 

It has been seen that there is no constant size difference between the 
gametangia in dioecious Mucors nor constant vegetative characters by 
means of which it is possible to distinguish the two sexes. The attempt 
to determine which is male and which female by imperfect sexual reactions 
with heterogamic hermaphrodites has been unsuccessful. 

A rather extensive series of qualitative and quantitative tests indicates 
that the (+) races of Mucors react biochemically in a similar manner to 
the females of higher plants and the (—) races like males.1° 1412 13 The evi- 
dence from these tests suggests that our (+) sex is female and our (—) sex 
male and is in line with the evidence from imperfect hybridization reactions 
with all except one species of the heterogamic hermaphrodites. Until 
further evidence is available, it will not be possible to be certain regarding 
the homologies of the two sexes in Mucors and higher forms. Asa matter 
of convenience, however, we may designate the (+) and (—) races as 
tentatively female and male. 

The detailed data upon which the present summary is based are being 
published elsewhere. 


1 The investigation was carried on under joint support of the Carnegie Institution 
of Washington and Mrs. Walter S. Franklin, Mrs. Walter B. James and Walter Jennings, 
to whom the authors gratefully acknowledge their indebtedness. 

2 Burgeff, H., Flora, 108, 353-448, 1915. 

3 Blakeslee, A. F., Cartledge, J. L., Welch, D. S., and Bergner, A. D., Bot. Gaz., 84, 
27-50, 1927. 

4 Blakeslee, A. F., Cartledge, J. L., and Welch, D. S., Jbid., 72, 185-219, 1921. 

5 Hartmann, M., Biol. Zentralbl., 45, 449-467, 1925. 

6 Blakeslee, A. F., and Cartledge, J. L., Bot. Gaz., 84, 51-57, 1927. 

7 Blakeslee, A. F., Biol. Bull., 25, 87-102, 1915. 
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§ Wilson, E. B., The Cell in Development and Heredity, 3rd ed., Macmillan Co., 
p. 592, 1925. 

® Hartmann, M., Verteilung, Bestimmung und Vererbung des Geschlecht ‘bei den 
Protisten und Thallophyten, Handbuch. Vererbung., B2, Lief 9, Gebriider Borntraeger, 
Chap. 2, p. 74, 1929. 

10 Satina, S., and Blakeslee, A. F., These PRocEEDINGS, 11, 528-534, 1925. 

11 Satina, S., and Blakeslee, A. F., Zbid., 12, 191-202, 1926. 

12 Satina, S., and Blakeslee, A. F., Jbid., 13, 115-122, 1927. 

13 Satina, S., and Blakeslee, A. F., Jbid., 14, 308-316, 1928. 


INFLUENCE OF ALUMINIUM ON MORTAR STRENGTH 
By H. W. Leavitt, J. W. GOWEN AND L. C. JENNESS 


MAINE TECHNOLOGY EXPERIMENT STATION, No. 8 


Communicated July 19, 1929 


The writers have for some years been analyzing the variation in the 
strength developed by mortars made from native sands and common 
commercial cements. This variation in actual test condition lies between 
100 and 600 pounds per square inch in tension and 2500 to 9000 pounds 
in compression. Among the constituents shown to account for part of 
this variation in strength is the amount of native iron* found in the sand, 
a correlation of 0.38 + 0.03 existing between the content of this native 
iron and strength. In this paper another sand constituent, aluminium, 
is studied and its relation to the strength developed by the mortar con- 
taining it. 

Table 1 shows that there is a consistent increase in the tensile strength 
developed by mortar on the 28-day test as the content of aluminium 
found in the native sand increases. This increase amounts to more than 
100 pounds over the whole range of variation in the native sands’ aluminium 


content. 
TABLE 1 


AVERAGE 28-Day TENSILE STRENGTH OF MORTARS IN RELATION TO THE VARIATION 
IN ALUMINIUM CONTENT OF THE NATIVE SANDS 


ALUMINIUM CONTENT, MEAN MORTAR STRENGTH, 
% NUMBER OF SANDS 1: 3 Mix 
0.25-0.50 17 298 
0.50-0.75 28 339 
0.75-1 .00 15 333 
1.00-1 .25 17 379 
1.25-1.50 15 385 
1.50-1.75 11 399 
1.75-3.75 14 418 


The amount of variation which exists for any one class and the degree 
of correlation which exists between the aluminium content and strength 
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may be seen in tables 2 and 3. 
Table 2 represents the 28-day 
tension results and table 3 the 
7-day tension results. 

Tables 2 and 3 show that as 
the amount of aluminium in 
the sand increases the break- 
ing strength increases. The 
straight lines on the tables 
show the general trend in the 
relationship. The correlation 
coefficients for aluminium and 
strength are: 7-day, 0.36 + 
0.05, and 28-day, 0.41 + 0.05. 
The correlation is, therefore, 
significant. It will be further 
noticed that within any given 
class of aluminium per cents, 
breaking strengths still have a 
fairly wide variation. Thus 
while the aluminium content 
significantly affects strength it 
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TABLE NO. 2 
Relation between aluminium content and 28-day 
strength. 


is clearly not the only limiting variable. In fact, we already know another 
chemical constituent which also plays a part, namely, iron. The part 


°7e Aluminium 


TABLE NO. 3 
Relation between aluminium content and 7-day 
tensile strength. 





which this iron content plays 
i in relation to the aluminium 
content will be analyzed in the 
/| succeeding paper. 

The amount of variation in 
the strength for which the vari- 
ation in aluminium content is 
responsible is shown by the 
correlation coefficients to be 
about 8% of the total variation 
in mortar strength. The fact 
that the correlation coefficients 
are approximately the same for 
the 7-day and 28-day tests 
shows that the aluminium con- 
tent continues to be equally ef- 
fective in raising or lowering the 
tensile strength throughout this 
duration of life of the mortar. 
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Aluminium content does not have the same marked effect on the com- 
pression strength of these mortars. In fact, the correlation between 
aluminium content and compression strength for these 117 sands is only 
0.115 + 0.062, scarcely significant. Thus again we have called to our 
attention a fact emphasized elsewhere in our papers, i.e., that the tensile 
test for mortar strength measures different attributes of the mortar 
strength than the compression strength. 

Summary.—The evidence presented clearly indicates the importance 
of another chemical element of native sand in affecting mortar strength. 
Aluminium in these native sands affects 7-day and 28-day tensile strength 
in such a manner that the larger its amount the greater is the tensile 
strength. Aluminium does not materially affect the strength of the test 
cylinders in compression. 

* “Influence of Iron Content on Mortar Strength,” by H. Walter Leavitt and John 
W. Gowen, Proc. Nat. Acad. Sci., 13, No. 4, pp. 263-265, April, 1927. 


ON THE JOINT INFLUENCE OF IRON AND ALUMINIUM IN 
NATIVE SANDS ON MORTAR STRENGTH 


By H. W. Leavitt, J. W. GOWEN AND L. C. JENNESS 
MaIn TECHNOLOGY EXPERIMENT STATION, No. 9 


Communicated July 19, 1929 


In the preceding papers the influence of iron and aluminium content 
of native sands on the strength of these sands in mortar has been analyzed 
separately. In this paper the joint influence of these elements on strength 
will be shown in its true perspective. Table 1 shows the relation of the 
iron content to the aluminium content of these 117 sands. Measured in 
terms of correlation this relation is 0.384 + 0.054. The effect of the 
iron content on mortar strength has been indicated elsewhere! to be 
0.377 + 0.034 for the 28-day period. A rapid digestion technique was 
used for the determination of the iron in the 284 sands involved. For 
the 117 sands here utilized the relation was 0.323 + 0.055 where the iron 
content was determined by a digestion test of 5 hours in hot HCl. The 
results are essentially the same. The aluminium content in relation to 
strength is shown in the preceding paper to be 0.412 + 0.052. 

From these three possible correlations the influence of any one of the 
elements on another with the third constant may be found. Thus the 
correlation for aluminium and strength with the amount of iron constant 
is 0.330 + 0.055, for aluminium and iron for strength constant is 0.291 + 
0.051, and for iron and strength for aluminium constant is 0.196 + 0.058. 
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The most important of these elements in determining the strength of the 
117 native sands is the aluminium content. Iron, however, still retains 
an important influence even with the aluminium constant. 

The joint influence of the two elements, iron and aluminium, is measured 
by the total correlation coefficient, r. This total correlation is equal to 
+0.538. The joint influence of iron and aluminium thus accounts for 
about 16% of the variation in mortar strength. 


form, the equation for the prob- 
able strength of the mortar deter- 
mined from the aluminium and the 
iron. This equation is: 


Mortar strength = 20 Aluminium% 
+ 43 Iron % + 272. 


The equation for strength shows 
that the two variables both in con- 
junction with each other and sepa- 
rately play a part in the strength of 
the mortar. The results therefore 
point to the conclusion that differ- 
ent chemical constituents of the 
native sands, aluminium and iron, 
affect tensile strength in mortars. 

Summary.—The results herein 
presented furnish proof of the fact 
that iron and aluminium jointly 


the sand and the cement. 








The influence of the two elements may be expressed in slightly different 
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TABLE NO. 1 


Relation between the content of iron and 
aluminium in native sands. 


and separately materially influence the strength developed by mortars 
made from native Maine sands. This conclusion we feel justified in 
extending to the more general type, the chemical condition of the sand 
used in cement mortars materially influences the strength of these mortars 
when hardened. Proof is thus furnished of the chemical interaction of 


1 “Influence of Iron Content on Mortar Strength,’’ by H. Walter Leavitt and John 
W. Gowen, Proc. Nat. Acad. Sci., 13, No. 4, pp. 263-265, April, 1927. 
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PLANETARY MOTION IN A RETARDED NEWTONIAN 
POTENTIAL FIELD 


By Roy J. KENNEDY* 
DEPARTMENT OF Puysics, CALIFORNIA INSTITUTE OF TECHNOLOGY 


Communicated July 31, 1929 


The problem of the motion of a planet within a solar system which is 
itself in uniform motion with respect to an ‘“‘absolute’’ space has been 
discussed from the time of Laplace up to shortly after the advent of 
Einstein’s gravitational theory. The latest work seems to have been 
done by Lodge! and Eddington,' who conclude that the pre-relativity 
discrepancies between theory and observation are not explicable in terms 
of variability of mass with velocity of the planet. 

Apparently no one has thought of applying to gravitation the idea 
corresponding to that of effective charge which Lienard? and Wiechert’ 
developed in electromagnetism. ‘The application is simple; in order to 
determine the potential of a moving gravitating mass we proceed as 
follows: 

The potential at any point 0 is given by the expression 


o=n ff eas 
le 


integrated over the boundaries of the body, dV, and r, being the ‘‘effective’’ 
volume and radius; the first is the element of volume in a fixed reference 
system within which each point of the moving element was located at 
the time when it gave rise to the potential existing at 0 at the slightly 
later time when the distance from the moving element to 0 is 7, while r, 
is the distance of the center of dV, from 0. Hence in rectangular fixed 
axes with x-axis in the direction of motion of the mass, if x, y, z, be the 
instantaneous coérdinates of a point of any element of the mass, and 
Xey Vey Ze the effective codrdinates, then i 
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where v is the velocity of the body and # is the time of propagation of 


gravitation from the point to 0. If the coérdinates of 0 are x’, y’, 2’, 


r = ct? (where c is velocity of gravitation) 


(x! — x)? + (9 — 9)? + (2’ — %)? 
(x’ — x)? + (y’ — vy)? + (2’ — 2)? + Qot(x’ — x) + vt. 


Or (c2 — v*)t? — Qv(x’ — x)t —- rr? = 0, 
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Xe = x — Br— B(x’ — x) (if 6% and higher powers are 
neglected) and 


ing « aw + Se de = (1+ Bcos6 + B2)dx 
Ox r Ox 


where @ is angle between x-axis and r,. 
. dV, = dy, dz, dx, = (1 + Bcos@ + B*)dV. 
Also = ct = r{1 + Bcos@ + 8? (1 + cos? 4)] 


eens p(1 + B cos 6 + B*)dV 
il 7 


1+ Beos0 + 5 (1 + cost 6) | 


wb ff f2(.+% sinto)av. 
r 2 


if powers of 8 above the second are ignored, as they will be henceforth 
without mention. 

Two distinct cases will be dealt with in determining the orbits; in the 
first it will be assumed that mass of planet is simply a constant independent 
of velocity and of direction of acceleration, while in the second the mass 
will be assumed to correspond to that of the Lorentz electron. In both 
cases the central mass can be treated as a homogeneous sphere, since only 
small corrections are being computed; hence, we may write for the present 


purpose 
= H+ + — Br sin? s) 
where yp is a constant. 


Case I.—Mass of planet constant. 
Here we can apply Lagrange’s equations in usual form 


d ar) oT 
£f ie oe ee: 1.1 
dt (3. . We Q (1.1) 


Choosing polar coérdinates 7, @ in plane of orbit, and using 6 for ratio 
of component of ‘‘absolute’’ velocity of the sun in that plane to velocity 





of gravitation, we have, since T = 5 (r? + 7176?) 
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tt = mr? ; or = 0 
of = mr : oT = mo? 
Q = ~ = 1 sin 20 (1.2) 
a. = — £(1 re sin’ ) (1.3) 


Substituting these values in (1.1) we get the simultaneous equations 


ad 4, _ »f* sin 20 

di 19) = > ae ge and (1.4) 
oo. ay. -# ee 
7) r (*) gaia (1 + > sin ) (1.5) 


Multiplying (1.4) by 278 we get 
@ (r26)? = pB’r sin 20. 6. 
dt : 
Now, approximately, 
h? 1 


r=— (1.6) 
uw 1+ ecos (@ — w) 





where h and w are constants, and e is eccentricity of orbit. 


_@ (9)? = h? 6? sin 20 dé 


dt 1 + ecos (6 — w) dt 





whence, on integration, 


2 3 
@= ak (f MBloadinad ) = oA. say. (1.7) 
Se x 1 + ecos (@ — w) r? 





Hence the operator 


_ dr _ hid (5 ) _ hy Ee 1 didr_ 21 G ' 
r? 0? © 2rIt dodo = r® 


‘de rt dé\r? db}? d6 


Substituting right side of this equation in (1.5) we get, after rearranging, 
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dr , 1dIdr_2(dr\? sr Bt, 
doi + 3T 0 do 2 (%) Fes — 1 (1 + & sint 4), 


Changing the dependent variable in this equation by setting u = . so 
dr 1 du d*r 1 du 2 (du\? 
eH do ~ ut db aod Te - — ge t (Ge) 


we obtain, on reduction, 


2. 2 2 
eum (14S sinte - Oe) 
n 





de? h?I 2 2u dé dé 
By (1.7) 
ay B? sin 20 du i oe 
aise ee — 7a é sin (6 w), 


to a sufficient approximation. Hence the differential equation becomes 


d*u : 6B? sin 20 sin (6 — 2) 
— 1 29 . (1.8 
a mal a ph 1 + ecos (6 — w) et 





When it is recalled that for an undisturbed Newtonian orbit, 6 = h, 
it is evident from (1.7) that 


aoe 
r=it f° B? sin 20 dé 
1 + ecos (6 — w) 


Hence (1.8) may be written 





du 
do? 


Be sin 20 sin (@ — w) 
2 1+ ecos (0 — w) 


~ p y sin 20d6 | (1.9) 
» l+ecos(@—w)} -— 
Solution of an equation of this form, i.e., 


ad? 
at 8 I 





+u m1 +E cinta 4% 





is 
u=Acos#+Bsino+ f° f(é sin (0 — &)dé. (1.10) 


In expanded form, f(@) evidently consists solely of terms in sin"@ cos”6 
so that solution will be derived from terms 


Ssin"é cos™tsin (0 — §)dé = sin 0f” sin" cos” *'tdé 
— cos 0f” sin" t' cos™édé. 
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Such integrals yield purely periodic terms (which contribute no observ- 
able effects) unless one of the exponents » or m is even while the other 
is odd. Furthermore, terms containing powers of the eccentricity e 
above the second are negligibly small. Thus the significant terms on 
the right side of (1.9) are 


2 
7 sin 20 sin (@ — w) + Be f-’sin 20 cos (8 — w)déd = f (8). 


Using this result in (1.10) and discarding the small periodic terms as they 
appear, we get as equation of orbit, 


B7e6 
u = B[ 1 + ecos 0 ~ «) — 2 sin (@ — 0) | 
2 
= "11+ ecos (@ — w) a F* og 20:0 0in (@ — w) 
h? 8 
2 
_ Fe sin 2w.6 cos (9 — »)| 


which may be written with sufficient approximation 


- F sin 20.8) (+- aes | 
u = Bl + (¢ 8 sin 2w.@] cos | 0 w+ gq COs 2.8 . 


Hence the perihelion apse advances at rate 


9 
bwo = — md B cos 2w 
8 
: : WAS ae ee 
per revolution, while eccentricity increases by 6@ = — 3 B’e sin 2w per 
revolution. Evidently 5e9 = tan 2w.edwo. Hence, using values in 
‘ 8.8 
table‘ on right below for Mercury, tan 2w = — a 0.104. 
*. 20 = —6° or +174°. Since dw is positive, cos 2w must be nega- 
‘ ss ‘ Qn 1 
tive; .. 2vW = 174°. The advance per century dw = — 7 B? cos 2w X = 


(where 7 is period of revolution expressed in years), or edw = 
257eB? cos 2w 


F Hence for Mercury, for which e = 0.206 and T = 0.24, 


3.48 2 3 ie 
gt =m — 848 X02 X5 X10" | ops + 10. 





25x X 0.206 X cos 174° 
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.. B = 7.97 X 10~‘ and if we assume the velocity of gravitation to 
equal that of light, the velocity of system is v = Bc = 7.97 XK 10-4 X 3 
xX 10+> = 239 km. per sec. This is by no means unreasonably large. 

Computing edw and de for the other inner planets on basis of above 
values for Mercury and table of longitudes of perihelion, and subtracting 
results from outstanding pre-relativity values of discrepancies we find 
as remaining discrepancies from observations per century: 


PRESENT THEORY NEWTONIAN THEORY 

€ bw be edw be 
Mercury 0.00 + 0.48 0.00 + 0.50 8.48 0.438 —0.88 + 0.50 
Venus —0.05 + 0.25 0.09 = 0.31 -0.05 +0.25 +0.21 + 0.31 
Earth —0.038 = 0.13 —0.10 = 0.10 0.10 + 0.13 0.02 + 0.10 
Mars —1.18 + 0.34 0.05 + 0.27 0.75 + 0.34 0.29 + 0.27 


Thus the present theory is in good agreement with observation, the 
discrepancy in perihelial motion of Mars being the only one greater than 
the probable error. Possibly the application of this theory to the per- 
turbation of Mars by massive Jupiter would remove the latter discrepancy. 

Case IIT.—Mass of planet variable. 

Here it is assumed that the inertial mass of the planet is dependent 
on direction of acceleration in same way as in the case of the Lorentz 


electron, i.e., the mass for acceleration in direction of motion is a 
while for acceleration in perpendicular direction it is a where 8 
is ratio of absolute velocity of planet in orbital plane to velocity of gravi- 
tation; the gravitational force depends only on the constant mo. ‘The 
treatment is a little less simple in this case, since the usual form of La- 
grange’s equations is not applicable. 

If rectangular axes are chosen with sun at origin, abscissa in direction 
of projection of absolute constant velocity vo of sun, and xy-plane in 
plane of orbit, and if 6, y, w and are angles from x-axis to radius vector 
to planet, direction of force, direction of perihelial apse and direction of 
instantaneous absolute motion of planet; respectively, then 


moo = f,(1 — 6’)? and moa, = fy (1 — 6)”, 


where a and f refer to acceleration and gravitational force, and subscripts 
to directions £ and perpendicular to &. 


T 
Mot, = Mo E cos~ + a, Cos (é + =) 
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Also fr = feos (§—g) and 
fn = —f sin ( — 9): 
Hence 


moa, = f{(1 — 6%)" cos = cos (— — y) + (1 — 6%)” sin & sin (E — ¢)] 


=H{(1 -£) [cos — cos (§ — ¢) + sin — sin (E — ¢)] 
— B? cos — cos (E — a} 


= if (i _ 2B Neos g — B’ cos ~ cos (E — | 
Similarly, 
B? 
May = (1 _ e) sin g — 6% sin — cos ( — | 


moa = m(o2 + a2)'/? = ak ~ us (1 + 2 cos*(é — o} 


Also 
a, _ cos o(1 — 36?) — 6? cos é cos (— — ¢) 


a sin g(l1 — 36?) — B* sin £ cos (§ — ¢) 





A 
cot xa 


cos (& — 


= cot » + 6? ¢) (sin £ cot g — cos &)- 


sin 9 


Since directions of force and acceleration are but little different, we have 


“ : Av; ; 
approximately on expansion of cot xa in Taylor’s series 


A ’ d 
cot xa = cot » + (tx — 9) “cot et... 
y 


A 
cot g — (xa — gy) csc*g +... 


ll 


Comparing the last two equations we find for angle from force to accelera- 
tion vector 


A 2 _ 
pe ae boa _ St (sin — cot g — cos &) 
sin ¢ csc? ¢ 
Bg? B? ; 
i sin 2(§ — yg) = — - sin 2 (§ — 0), approximately. 
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Also it follows from (1.2) and (1.3) that angle from force vector to radius 
vector is 


> 
6@-—g= 5 sin 20 
where By = =. Hence angle from radius vector to acceleration vector is 
A af Bo 
v = xa —60= (xa — ¢) — (0—¢) = 5 sin 20 — £) — 3 sin 26. (2.1) 


Now since corrections to orbit are very small we can of course compute 
8 from velocity in a truly elliptic orbit. By an elementary theorem the 
orbital velocity may then be resolved into a constant component ev 


inclined at constant angle w + 5 with x-axis and another constant com- 


ponent v, perpendicular to radius vector, i.e., inclined at angle @ + 5 with 


x-axis. Hence the orbital component of absolute velocity v,’ of planet 
has components 


Vx’ = Up + ev, COS (. + 5) + v cos (6 “he *) 


= UW — ev, Sin w — 2 sin 0 and 
Vy = ed; COS w + % Cos 4: 


id , 
2 
v? uw +r 


c c? 


. B? 
= B+ 5 { [vf(1 + e?) — 2ev,v sin w] + ev} cos (0 — w) — 2% sin O}. (2.2) 


(2.1) may be written 


2 is 2 
bel: [sin 26(1 — 2 sin? £) — 2 cos 26 sin — cos —] — fs sin 20. (2.3) 


Now 
, , 
vy 0 
sin = ~ = — (cos 6 + e cos w) 
v cB 


and 


cos @ = = "11 — "(sino + ¢ sin w) | 
v cB 


Vo 
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Substituting these values in (2.3), and then for 6? from (2.2) we obtain 
after a few transformations 


y= — a [cos 8 + e cos (20 — w)] 





Vv 


2 
ae le sin (0 — w) + < sin 2(6 — »)| (2.4) 
c? 2 


If component of acceleration perpendicular to radius in direction of 
motion is ag, we have 


T : e 
a = 200s (7 +7) = —asinvy = —a=—- v= eee (255) 


to usual approximation. 
Similarly the acceleration along radius is 


a, = acosyv = -An Fir t 2eost@ ait - 7) 


v is negligible. Expanding cos* (@ —£) and substituting above values 
of sin &, cos and 6? and dropping small constant terms and terms in 
sin” @ cos” 6 for which both m and n are even or odd, we find 


U1V0 
C2 





a = — ar + sin 6 — “te cos (6 — w) | (2.6) 


In general, the components of acceleration along and perpendicular 
to radius vector are, respectively, ry — 76? and ré + 276. Equating these 
values to (2.6) and (2.5), 


in 0 6 cos (0 — | 


C2 





2 


r— 162 = -Ali+ 
n 


S 


and 


r6 + 279 = — Hy. 


_ 


Treating these equations just as the corresponding ones in Case I were 
treated, we obtain the differential equation 


hi? dl du 
2u dé dé |’ 


d*u i U0. V} 
Te um [1+ sino — Le cosa — w) — 
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=f — vdé -1-[ vd6 
1 + ecos (@ — w) o 1+ ecos (6 — w) 


where 








Hence finally, if terms containing third and higher powers of e are dropped, 


010 ; 


Hunt efi + Mss sin @ — “Se cos (0 — w) 


e vsin (@ — w) ft ; vd@ 
— = + 
21+ ecos (0 — w) o 1+ ecos (@ — w) 








Solving this equation by method used with Case I, we obtain 


u = Bt + e008 0 — 0) +5 2 005 0.8 cos @ — «) 
+ 3 of sin w.0 sin (0 — w) +5 Me cos w.6 
16 c* 
~4§ Se pu ems 
Bhit(e+5e a cos w.8) cos (# w 16° c sinw.) 
+ 5 —F cos w.0 


to approximation previously used. 

Thus the orbits will oe changes of eccentricity and perihelial phase 
per revolution of 27. = e? ae cos w and + 2r. a e — sin w, respectively, 
together with a slow variation in radius. If w is computed from 
the observed motion of Mercury, it is found that the elements for other 
planets are nearly zero; the agreement with observation is therefore good. 
The rate of variation of radius even for Mercury is small, somewhat below 
the limit of observation. However, the value of velocity of solar system in 
its own plane turns out so large (about 18,000 km. per sec.) that it would 
surely have shown itself in discrepancies in the times of transit of Jupiter’s 
moons; the value of 7% computed in Case I is just too small to be observ- 
able in this way. 

In view of the requirement of covariance set forth by Bateman and 
Einstein, a theory of the kind developed in this paper can scarcely be 
given great weight; nevertheless it does seem to show that the classical 
poteritial theory is capable of accounting for planetary motion about as 
accurately as does the relativity theory, if less probably. 

* Guggenheim Fellow in Physics. 
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ON THE VIBRATIONAL SELECTION PRINCIPLES IN THE 
RAMAN EFFECT 


By J. H. VAN VLECK 
DEPARTMENT OF Puysics, UNIVERSITY OF WISCONSIN 


Communicated August 13, 1929 


Various writers’®* have noted that the displacements of energy levels 
involved in the Raman effect should be those obtainable by superposition 
of two ‘‘allowed” spectral transitions. In other words, v(ac) is a possible 
Raman shift if we can write v(ac) = v(ab) + v(bc), where v(ab), v(bc) are 
possible emission or absorption frequencies. The reason for this is that 
the theoretical expression for a Raman line of frequency, v + v(ac), excited . 
by incident light of frequency » is proportional to 


(vy + v(ac))*{| XP ee) |? + [YP eg |? + |Z ee) |*}, (1) 


E yp, [Plabelbe) _ ob PY) patna, 
vo + v(bc) vo + v(ab) 


Eqs. (1, 2) are the same as those recently given by Hill and Kemble,’ and 
embody the results of the Kramers dispersion theory. Here X(ab).. 
denote the matrix elements of the Cartesian components of the electrical 
moment of the molecule in the unperturbed state, while x® (ab)... . denote 
the perturbed or ‘““Raman’”’ values of these elements under excitation by 
plane polarized incident radiation having its electric vector E cos 2rv ¢ 
along the z direction. Eq. (2) involves the amplitudes connected with 
transitions to what we shall term the ‘intermediate’ states b rather than 
the amplitudes P(ac) of direct transition from a to c, thus substantiating 
the remarks made at the beginning of the paragraph. 

Notation, etc—Before even stating the particular problem connected 
with Eqs. (1, 2) which we aim to discuss, it is perhaps well to introduce 
some nomenclature and qualifying remarks. We shall suppose except 
in the final paragraphs that the molecule is diatomic, and shall ignore 
the rotational fine structure. The letters v and n will denote, respectively, 
the vibrational quantum number and the totality of electronic quantum 
numbers. To separate the vibrational and electronic effects, it is often 
convenient to use a double index n, v for the matrix arguments instead of 
a single index such as a, b, c in Eq. (2). The term “‘level” is to be under- 
stood to comprise all the vibrational states having a given set of electronic 
quantum numbers, and we shall use the term ‘‘state’’ for the individual 
vibrational members of a level. The normal level will be denoted by m. 
Two levels a, b will be said to “‘combine’’ if the unperturbed amplitude 
P(ab) does not vanish; i.e., if v(ab) is a possible absorption (or emission) 





where 





(R) as 
P (ac) ~ 
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frequency, but not necessarily a possible Raman shift. We shall use the 
term “changes” in discussing the alterations of quantum numbers in 
ordinary absorption or emission transitions, and shall reserve the word 
“displacements” for the total alterations in quantum numbers between 
the initial and final Raman states. It is to be clearly understood that the 
term ‘‘intermediate’’ relates merely to the position of a state such as b 
in the products in Eq. (2), and does not at all imply that its energy is 
intermediate between those of the initial and final Raman states, as this 
is usually not the case. 

Problem of the Small Raman Displacements in the Vibrational Quantum 
Number.—The most commonly studied Raman lines represent displace- 
ments in the vibrational, but not in the electronic quantum numbers. 
However, the important intermediate states involved therein probably 
belong to excited electronic levels, as the absolute intensities of pure vibra- 
tion spectra are low and hence the amplitudes connected with transitions 
to intermediate vibrational states belonging to the normal level are small.‘ 
Thus Raman displacements corresponding to infra-red vibrational fre- 
quencies are due in a certain sense to excitation of electronic absorption 
lines usually in the ultra-violet. The incident light is ordinarily in the 
visible region, and so could not excite these lines with close resonance, 
but nevertheless does react somewhat with the amplitudes of ultra-violet 
absorption, affording an illustration of what has sometimes been termed 
“weak quantization,” especially in the old quantum theory. Now when 
there are electron jumps, the selection principles for the vibrational 
quantum number are quite different from what they are in pure vibration 
spectra. Whereas we can only have Av = + 1 in simple harmonic, 
pure vibration motions, there can be enormous changes in v if there are 
simultaneous changes in n. In the latter event cases are known in which 
v changes by 20 units or so.» Now if a Raman displacement is expressible 
as the result of two consecutive jumps, one might at first thought expect 
Raman displacements of say 20 + 20 or 40 in the vibrational quantum 
number (taking this rather extrer:e example for concreteness). 


More precisely, the allowed changes in v are determined by the requirement that the 
initial and final states are to have the same inter-nuclear distances during appreciable 
time fractions of their respective motion: This is the rule advanced semi-empirically 
by Franck® and justified quantum-meche tically by Condon. Now the nuciei move 
slowest and spend relatively the greatest i me at the extremes 7,,;,, 7max Of their vibra- 
tional paths for a given state. Hence two states combine intensely if they have a 
common extremum of the inter-nuclear distance 7. This extremum can, however, be 
a maximum for one state and a minimum for the other. In particular, if we start 
with a given vibrational state v of a given electronic level m, changes are particularly 
favored to two vibrational states of another electronic level n’ (viz. the Pair fain = main 
ry = Tmax if min < r bs Tmax) the pair Fane = Tmins foe = Tmax if r< Tmin the pair 


max 
ri, i ue ro >fmax; here ry denotes the equilibrium value of r in the 


wei: min 


min? 
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state n’). This, as Condon shows, is in good agreement with experiment. If one 
plots initial and final vibrational quantum numbers as abscissa and ordinates, the 
intense changes form a ‘‘Condon’”’ parabola, and the ordinate or abscissa representing 
a given state intersects this parabola twice, giving the two favored transitions from this 
state. To apply these ideas to the Raman effect let us start in Eq. (2) with a given 
initial level a = mo, vo belonging to the normal electronic level. Let us further consider 
the effect of all the intermediate vibrational states belonging to a given electronic level 
n. In view of the foregoing there will be two of these states, say m, v; and n, ve which 
combine particularly with mo, vo. If the Raman displacement is purely vibrational, 
the final state c belongs to the normal level m. But the intermediate state n, 1, can 
combine with another vibrational state, say mo, v3, of the normal level besides m0, vo, 
as there are two favored transitions. Similarly mo, ve can combine with another state 
No, v4 besides mo, vo. As a Raman displacement is compounded from two ordinary 
changes, one might expect intense Raman lines corresponding to displacements vo — v3 
and vp — vin the vibrational quantum number, besides of course the Rayleigh scattering 
Av = 0. Such displacements represent passage from one side to the other of the Condon 
parabola for the intermediate level. When finally one includes all the different inter- 
mediate electronic levels in (2) one would obtain an enormous range of possible dis- 
placements in v, of large amounts if the important intermediate levels have widely 
different moments of inertia from the level mo. 

These hasty conclusions are not in agreement with experiment. Actually 
the intense Raman lines, at least in gases, involve shifts in frequency equal 
to the fundamental line of a pure vibration spectrum, or perhaps in some 
cases the first harmonic. This means that the common Raman displace- 
ments in the vibrational quantum number are only one unit, or possibly 
sometimes two units. The flaw in the foregoing argument predicting large 
displacements is the following. One must remember that according to the 
Condon-Franck principle, a given state combines, though with greatly 
varying intensity, with all the vibrational states belonging to another 
electronic level. In particular, we assumed in the preceding fine print 
that a given state combined practically entirely with only a pair of states 
of another level, instead of with all other vibrational states of this other 
level, though with reduced intensity. The reduction in intensity is, as 
a matter of fact, not large for states which are immediate neighbors of 
either of the states (the favored pair) that combine most strongly with 
the given initial state; this is especially true if the vibrational quantum 
numbers are large. In forming the sum in (2), one must therefore include 
the contribution of every vibrational state belonging to a given inter- 
mediate level. Furthermore, in adding these contributions, one must 
be very careful to consider the phases of the matrix elements involved 
therein. One finds that the phase effects give reinforcement when the 
initial and final levels a and c are the same, the case of Rayleigh scattering. 
(At least this is true for the z component, as /(ab)/(ba) = | /(ab) |2 >0 
and so all the terms to be added have the same sign, as with a given inter- 
mediate level the variation in the denominator from term to term in the 
summation in (2) is negligible.) On the other hand, because of the di- 
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versity of phase factors for the different intermediate states belonging to 
a given electronic level, one suspects that the various terms in the sum- 
mation will largely annul each other if the initial and final states a and c 
differ widely in their value of v. If, however, the states a and ¢ are ad- 
jacent, i.e., if the Raman displacement in v is only one unit, the annul- 
ment might be very incomplete, as this case is very close to that of re- 
inforcement. By a somewhat similar qualitative reasoning, Hill and 
Kemble* briefly stated that the phase relations explain why there are 
only small Raman displacements in the vibrational quantum number. 
The purpose of the present paper is to give a more quantitative and 
rigorous proof of this proposition of the usual absence of all displacements 
except the fundamental Av = +1. The writer is informed that C. 
Manneback is also publishing a proof, but by a method somewhat different 
from ours, so that the present note is perhaps not superfluous.’ 

Mathematical Preliminaries—We shall first develop a multiplication 
scheme for what we shall call “hybrid matrix elements,’’ whose use is the 
essence of the proof. Let WV, and %, (v = 0, 1, . . .) be two diflerent sets 
of orthogonal normalized functions of a variable r, so that 


SVVy dr = 8, Sf 8,%,dr = 8, (3) 


with the asterisk denoting the complex conjugate. For example, and more 
especially for our purposes, VY, and ®, might be the vibrational wave 
functions for two different electronic levels. Let f(r) be any function 
of r, and also if desired of the operator 0/0r. Expand f(r), as a series 
in the ®,, thus 

f(r)¥, = Dv fow (v’; v)®,. (4) 


The coefficients in this expansion we shall term the hybrid matrix elements 
of f, as they differ from ordinary matrix elements in that the wave function 
on the left belongs to a different set of-orthogonal functions from those on 
the right. The ordinary elements correspond to the special case VW = ®. 
On multiplication of (4) by ®,- and integration throughout the domain 
of 7, we see, using (3), that 


fay(v"; 0) =f Oy f¥edr. (5) 


Let 9,(r) be still a third set of normalized orthogonal functions. The 
important property of the hybrid elements is that they obey the following 
law of multiplication 


(feow(v"; 1) = Div foa(v"; v’)gey(v’; v). (6) 


To prove (6), note that fg¥,(7) can be expressed equally well as 
Livr( fedew(v"; v)®,-(7), (7) 
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or as 


fdvSev &(v’; v) = sav »” foo(2"; v’)geu(v’; v)O,-(r). (8) 


The result (6) follows immediately on comparing coefficients in (7) and 
(8). Besides the above multiplication law, the hybrid matrix elements 
obviously obey the ordinary elementary addition law. Incidentally, these 
elements for f = 1 are not the ordinary Christoffel symbols 8 since 
and ® are not orthogonal, but do have the property that 


Dv'les(v"; 0’) leg (v’; v) = <. 


As we have agreed to disregard the rotational factor, the total wave 
function of the molecule is approximately the product of an electronic 
wave function U(x, . . . 2;,7) and a vibrational one R(r). Here R is a 
function only of the inter-nuclear distance 7, while U is a function of the 
Cartesian coérdinates of the electrons, measured in a system of reference 
fixed relative to the nuclei, and also involves r as a parameter. As the 
electronic and nuclear variables are not rigorously separable in the Schroed- 
inger equation, this resolution of the wave function into electronic and 
vibrational factors is not exact, but is a good approximation if the vibra- 
tional distortion is not too large. It permits us to resolve the total energy 
W,,. into electronic and vibrational parts, so that 


Wr = W(n) + W,(2), (9) 


where W(n) is the energy for a system in which the nuclei are at rest 
at the equilibrium distance appropriate to the state n, and W,,(v) is the 
vibrational energy measured relative to an origin at this equilibrium 
position. By the Schroedinger-Eckart rule, the elements of any function 
f are given by*® 


f(n'v'; mn) =f... S [Un Raw f Un Ruy dor... de¢dr. (10) 


If we omit the vibrational wave factors and integrate over the electronic 
coérdinates, we thereby define a function of 7 


farnlt) =S 0. [ULF Ug dey... dey (11) 


which of course involves the initial and final electronic quantum numbers 
n', n as parameters. We shall write in the electronic quantum numbers 
as subscripts whenever they appear as parameters in the vibrational 
problem. Egs. (5, 10, 11) show that a complete matrix element’ is a 
hybrid matrix element in which f,,(7), Ruy, Ray now play the same roles 
as f(r),®,, ¥, in the notation used in connection with Eqs. (3-5). The 
rule (6) shows that 


(farm Bn’n)(0"; 0) = Dir Surnr (0; O’)Bnin(0"; 0), (12) 
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which is not to be confused with the ordinary law of complete matrix 
multiplication in which the summation over the intermediate states is 
with respect to both m and v, and so includes all possible electronic levels, 
whereas in (12) it includes only the vibrational states of one particular 
electronic level. It is to be understood that for this reason an element 
(v”, v) Of fun’ Sn’n iS not the same as the corresponding element of ( fg) n*n. 

Proof That Only Small Raman Displacements in the Vibrational Quantum 
Number Are Possible—After the mathematical preliminaries, this proof is 
quite simple. In view of (9) and the significance of the states a, b, ¢ in 
the more explicit notation, the Bohr frequency condition may be written 
as follows: 


hv(ab) me W (no) + Wro(V0) ree W(n) si W,(v'), 
hv(bc) = W(n) + W,(0’) — W(no) — W,,(v"). (18) 


Let us substitute (13) in (2), and for simplicity sum until further notice 
over only the intermediate states belonging to one particular electronic 
level. As the vibrational energy is small compared to the denominators 
in (2), we may expand (2) as a power series in W,,(v’), thereby bringing 
W,,(v’) in the numerator rather than denominator. Similar expansions 
in W,,(v) and W,,(v") prove unnecessary, as they are constant with 
respect to the summation. We thus find for the part of (2) resulting 
from the intermediate level n 


aDy {IA} — A Wy lo’) + A Wale)? +.) [Plorn; no") Z (0; no0")] 
— [B+ B? W,(v') + B*W,(0’)? +. .]Z | (10; nv’) P(nv; nv"), 
where (14) 


as Sat eer, A= hyo + W(n) arnt W (m) a W,,.(v”), 
B= hy + W(m) + W,,(v) — W(n). 


We have included only terms of the first and second orders in W,,(v’), but 
incorporation of higher orders would occasion no difficulty and would not 
affect our conclusions. Now the W,,(v’) are the diagonal matrix elements . 
of the Hamiltonian function H,,, for the vibrational problem of the state 
n, while the non-diagonal elements are zero as the energy is a diagonal 
matrix. Hence the first part of (14) can be written 


ardor’, gill { Pron (2; v’) (A oF ee | a + ao Hin) (0'; eee? ial v"). 
(15) 
Here we have used the notation P,,,,,(v; v’) rather than P(nv; nv’), etc., in 
order to conform to the usage in Eq. (12). We henceforth do not bother 
to write out the terms arising from the second part of (14), as the manipula- 
tion of these terms is similar to that of those from the first part, and can 
be supplied by the reader. The hybrid multiplication law (6) or (12) 
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when extended to a triple rather than double product, and use of formula 
(5) shows that (15) is the same as 


aS Row { Pron(A-? — A? Hun + A? H2y)Z nme} Rnoor dr. (16) 


Here the P,,,,, etc., are constructed in the fashion (11) by integrating only 
over the electronic coérdinates, and H,,,, is the operator 


1 h \? 0? 
ews a aa = a Ve ’ 
2M (5) an fa) 


corresponding to the vibrational problem of the state n. The important 
thing is that in (16) the summations over v’, v'! have been entirely elimi- 
nated. To evaluate (16) exactly it would be necessary to know the 
vibrational potential energy function V,,(r) for the excited or intermediate 
level n, but the vibrational wave functions K,,, for this level are no longer 
required. Only those Ry, Ry» of the normal level are needed. The 
dependence of the factors P,.,, Zn», on 7 is also required, and is in general’ 
not known in detail. We may, however, expand these expressions as 
power series in (r — 1) if the vibrational excitation in either of the states 
mv or mv" is small. For then either R,,. or Ry» is vanishingly small 
except in the vicinity of the equilibrium value ro of the inter-nuclear 
distance for the state m. This condition is presumably met in the initial 
state mv, as the concentration of molecules in excited vibrational states 
is usually negligible under ordinary conditions of temperature, etc. The 
expansion of (16) takes the form 


aS Rios { [a + b(r — m) + c(r — m)? +...][A-! — A-?H,, + A-3 H?,] 
[g + h(r — 10) + i(r — 70)2]} Raye dr. (17) 


Since the equilibrium distance 7, for the state m is not the same as that 
ro for the state m, the expansion of the potential energy for the state n 
will contain constant and linear terms in r — 1%, whereas the expansion 
of that for the state mp) will commence with quadratic terms. Hence the 
expansions of the vibrational Hamiltonian functions for these two states 
are, respectively, of the form 





1 


as ay oO et ae 
Hu = (+) —+d+elr ro) +f(r 1)", 


2ri/ Or? 
1 h \? 0? 
a = ee (+) 3 Y — 1)". 
taal ae * ~ 


Subtraction of these two expressions yields us a formula 


Haun 2 FL assy + d + e(r ae ro) + (f — fo(r si ro)? Piaiais (18) 
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for H,,,, which can be substituted in (17) and which is especially convenient 
because it does not involve the differential operator (h/2m7)?0? . . ./dr?. 
Unless some of the coefficients a, b, etc., happen fortuitously to be par- 
ticularly small, we can suppose that the terms in (17) having factors 
(r — ro)°, (r — 10), (r — 70)? are of rapidly decreasing orders of magnitude, 
as we supposed at least one of the wave functions Ry», Ryo» very small 
except in the vicinity of 79. To a sufficient approximation, the vibrational 
problem of the normal level may be supposed simple harmonic, making 
Riw Rn Hermitian orthogonal functions. Then by the familiar 
selection rule for the harmonic oscillator, the matrix elements of r — 17 
are entirely of the form Av = +1, while those of H,,,,, are, of course, 
entirely of the type Av = 0. (The elements of H,,,, would not be of this 
type, as they are diagonal with respect to the R,, rather than R,,,; for 
this reason the substitution of (18) in (17) isnecessary.) By the Schroed- 
inger-Eckart rule, the integral (17) is merely a matrix element of the 
part of the integrand inclosed in parentheses. Hence by the rules for 
matrix multiplication, the parts of (17) which are of orders (r — 19)°, 
(r — ro), (r — 70)? vanish unless, respectively, Av = 0, Av= +1, Av = 0, 


+2. Thus the most intense scattered lines are of the form Av = 0, 
i.e., represent simple Rayleigh scattering, while the most intense Raman 
lines are those Av = +1 whose displacement equals the fundamental 


of the vibrational frequency. Less intense still are the Raman lines 
Av = +2, whose displacement equals the harmonic of the vibrational 
frequency, and which are ordinarily far too weak to detect experimentally. 
This is the desired result. There can be no question of large vibrational 
transitions, corresponding to passage from one side of the Condon para- 
bola to the other for the intermediate level, as this idea has been dropped 
out of consideration by using the hybrid multiplication law. Our proof 
has hitherto involved the restriction of the summation in (2) to one par- 
ticular intermediate electronic level, but the inclusion of all levels®!® 
merely yields a linear combination of expressions of the type form (17), 
and so cannot cause large Raman displacements in v. 

Polyatomic Molecules.—Essentially the same considerations as in the 
rest of the paper can be applied to molecules composed of more than two 
atoms. The main difference is that instead of one vibrational coérdinate 
r, we must utilize several 7, 72, . . . which are identical with “normal co- 
ordinates”’ if the vibrations are simple harmonic, and which are associated 
with different kinds of vibrational quantum numbers 2, v2, ... The 
same general method as before is applicable. We find that the scattered 
lines should rank as follows in intensity: first, Rayleigh scattering; second, 
Raman lines whose displacements equal one of the fundamental vibration 
frequencies; third, Raman lines whose displacements are either harmonics 
of a fundamental vibration frequency, or else are “‘combination fre- 
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quencies” equal to the sum of two different fundamentals. However, 
Raman lines which are usually considered to be of the combination type 
have been reported experimentally without any trace of Raman shifts 
equal to the fundamental absorption frequencies. In other words, lines 
which we have ranked third in intensity are apparently found without 
those we ranked second. Instances in solids or liquids will not be cited, 
as our theory is primarily for gases. In CO, Rasetti? finds only Raman 
shifts of 1284 and 1392 (+ 10) cm.~! which are not equal to absorption 
lines, but which are expressible as difference combinations of the two 
components (3616 and 3720 cm.~') of the doublet absorption band at 
2.7u with that at 4.28u (2336 cm.—!). To account for this anomalous 
situation two different suppositions (I) or (II) might be made. 

(I) One possibility is that the intense Raman lines in CO, are funda- 
mentals, and have been misinterpreted as combinations. Eucken’s!'! col- 
linear model of the CO, molecule does indeed have a fundamental at 
1174 cm.~!, which Rasetti!? notes is in the vicinity of his Raman lines, 
but it is doubtful whether all the details of Eucken’s model are tenable. 
It is perfectly possible for a fundamental to be inactive as an absorption 
line, but to appear asa Raman line. This will be the case for symmetrical 
vibrations, which do not create an electrical moment in the molecule in 
the ordinary (dielectric constant) sense, but which will in general change the 
magnitude of the moment matrix elements associated with transitions to 
excited electron levels, and so cause non-vanishing coefficients b, h in 
Eq. (17). In Eucken’s model, for instance, the 1174 cm.~! vibration is 
inactive in absorption and corresponds to symmetrical displacements 
of the two O nuclei relative to the C nucleus at the center, thus merely 
stretching the molecule without altering its symmetry. The sole Raman 
displacement which Wood" finds in NH; vapor corresponds to an observed 
absorption band at 3u. The very careful recent analysis of the NH; 
vibration spectrum by Stinchcomb and Barker’ shows that the 3u vibra- 
tion is due to oscillations in the altitude of the molecular pyramid, which 
has.the N nucleus at its vertex and the H nuclei at the corners of its base. 
Their measurements do not reveal whether the 3u vibration is a funda- 
mental, but Wood’s observation of a 34 Raman displacement is strong 
evidence that this vibration is really a fundamental, and it is indeed so 
classed in Dennison’s'® model of the NH; molecule. The non-appearance 
of the other fundamentals as Raman lines suggests that stretching the 
altitude of the molecular pyramid distorts the moment matrix elements 
more than do the other normal modes of vibration. 

(II) In the particular’ case of CO, another possibility, for which the 
writer is indebted entirely to Professor Dennison, is that the Raman 
displacements observed by Rasetti are not pure vibration lines, but rather 
a combination of rotation and vibration: In a collinear model the energy 








—_ — = 
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of small oscillations of the C nucleus in the median plane (the perpendicular 
bisector of the line segment joining the two O nuclei) is proportional to 
2n, + n,.” Here n, and n, are respectively the radial and azimuthal 
quantum numbers for motions of the C nucleus in the median plane. 
Hence the apparent vibrational quantum number is 2m, + n,. Now 
Kemble and Hill,* for instance, have shown that only Raman displace- 
ments of only 0 and 2 units are possible in rotational quantum numbers 
such as ,. Hence in the Raman effect 2n, + n, can be displaced only 
by an even number of units. In investigations not yet published, Denni- 
son interprets the absorption band of CO, at 673 cm.~! as due to a unit 
change in 2n, + n,. The Raman displacements found by Rasetti are 
approximately double this, and hence might be interpreted as displace- 
ments of 2 units in 2m, + n,. This does not contradict the preceding 
theory in our paper because n, is a rotational rather than vibrational 
quantum number. It is, however, quite different from the usual rotational 
quantum number in that its effect on the energy is as large as that of vibra- 
tions. This is because of the anomalous collinear character of the CO, 
molecule, in virtue of which the moment of inertia about the axis of figure 
is very small, in fact, zero when the C nucleus is exactly at its equilibrium 
position. 
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